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We investigate how current and proposed observations of neutron stars can lead 
to an understanding of the state of their interiors and the key unknowns: the typical 
neutron star radius and the neutron star maximum mass. We consider observations 
made not only with photons, ranging from radio waves to X-rays, but also those 
involving neutrinos and gravity waves. We detail how precision determinations of 
structural properties would lead to significant restrictions on the poorly understood 
equation of state near and beyond the equilibrium density of nuclear matter. 

To begin, a theoretical analysis of neutron star structure, including general rela- 
tivistic limits to mass, compactness, and spin rates is made. A review is the made of 
recent observations such as pulsar timing (which leads to mass, spin period, glitch 
and moment of inertia estimates), optical and X-ray observations of cooling neutron 
stars (which lead to estimates of core temperatures and ages and inferences about 
the internal composition), and X-ray observations of accreting and bursting sources 
(which shed light on both the crustal properties and internal composition). Next, we 
discuss neutrino emission from proto-neutron stars and how neutrino observations of 
a supernova, from both current and planned detectors, might impact our knowledge 
of the interiors, mass and radii of neutron stars. We also explore the question of 
how superstrong magnetic fields could affect the equation of state and neutron star 
structure. This is followed by a look at binary mergers involving neutron stars and 
how the detection of gravity waves could unambiguously distinguish normal neutron 
stars from self-bound strange quark matter stars. 
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I. INTRODUCTION 

Neutron stars are the most compact known objects without event horizons and therefore serve as 
extraordinary laboratories for dense matter physics. The internal composition of the cores of neutron stars 
is currently poorly understood. Most models of dense matter predict that above twice the equilibrium 
density of nuclear matter, p s ~ 2.7x 10 14 g cm -3 or n s ~ 0.16 baryons fm -3 , exotica in the form of hyperons, 
a Bose condensate of pions or kaons, or deconfined quark matter, will eventually appear. However, whether 
the threshold density for such exotica is around twice p s , or much larger, is unclear. 

It has been suggested that if strange quark matter is the ultimate ground state of matter (i.e., has 
a lower energy at zero pressure than iron) [l|, compression of matter to sufficiently high density would 
trigger a phase transition converting virtually the entire star into strange quark matter. Such a star is 
self-bound as opposed to being gravitationaly bound as is the case of a normal neutron star. 

It has so far proved very difficult to find venues from astrophysical observations that could unambigu- 
ously distinguish strange quark stars from normal neutron stars. This is because self-bound stars have 
similar radii, moments of inertia, and neutrino emissivities and opacities to that of moderate mass normal 
neutron stars. Therefore, it may be unlikely that photon or neutrino observations, or radio binary pulsar 
timing measurements, will be able to differentiate these cases, especially if strange quark stars have a small 
hadronic crust, supported perhaps by electrostatic forces or a mixed phase 0. In that case, the effective 
temperatures and radii of solar-mass-sized strange quark stars and normal neutron stars would tend to 
be similar. Even during the proto-neutron star stage, which is observable through neutrino emissions 0], 
these two types of stellar configurations likely yield similar neutrino signals until such late times that the 
low luminosities prevent an unambiguous discrimination However, major differences in the evolutions 
of normal neutron stars and strange quark matter stars emerge during the final stages of binary mergers 
if stable mass transfer occurs. These differences would be prominent in both the amplitudes and frequen- 
cies of gravitational wave emisions. Except for discussing this latter topic, we will not discuss self-bound 
strange quark matter stars in detail further in this paper. 

Possibly the two most important properties of neutron stars - their maximum masses and typical radii 
- are not yet well known. These properties reflect rather different aspects of the dense matter equation 
of state (EOS). The neutron star maximum mass, which is a consequence of general relativity and does 
not exist in Newtonian gravity, has a limit of, at most, 3 M , assuming causality (5[. The maximum mass 
is controlled by the stiffness of the dense matter EOS at densities in excess of a few fimes times n s . The 
introduction of non-nucleonic degrees of freedom at supra-nuclear densities generally implies a softening 
of the EOS. The largest precisely known neutron star mass is only 1.44 M , which is too small to provide 
an effective constraint on the composition of dense matter and the structural properties neutron stars. 
However, some recent mass measurements, detailed below, from timing of pulsars in binaries with white 
dwarf companions suggest that the maximum mass might be about 2 M or even larger. In addition, large 
masses are also suggested from observations of QPO's, or quasi-periodic oscillations, observed in X-ray 
emission from neutron stars accreting from binary companions. If verified, a maximum mass in the vicinity 
of 2 M Q would fundamentally alter our perception of the properties of matter in neutron star cores. As 
we demonstrate below, sufficiently large observed neutron star masses set interesting upper limits to the 
densities possible in neutron stars. It could be the case that the maximum density allowed is not large 
enough for exotic forms of matter to play significant roles. 

On the other hand, the neutron star radius is controlled by properties of the nuclear force in the 
immediate vicinity of n s , in particular by the density dependence of the nuclear symmetry energy [6| (the 
symmetry energy is the difference, at a given density, between the baryon energies of pure neutron matter 
and symmetric nucleonic matter). Measurements of the neutron star radius are far less precise than mass 
measurements. Among the observations recently utilized include radius upper limits from rapidly rotating 
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neutron stars, estimates from the thermal emission of cooling neutron stars, including redshifts, estimates 
obtained from the properties of sources with bursts or thermonuclear explosions on their surfaces, and 
estimates of crustal properties from (a) glitches of pulsars, (b) "star-quakes" occurring in the aftermath of 
giant flares from soft-gamma ray repeaters (SGR's), and (c) cooling timescales during periods of quiescence 
in-between X-ray bursts from accreting neutron stars in low-mass X-ray binaries (LMXRB's). 

Aside from direct mass and radius determinations, one of our best windows into the interior is through 
observations of their thermal properties. Neutron star cooling curves (i.e., their luminosities or tempera- 
tures as a function of age) are sensitive to the internal composition and the superfluid characteristics of 
its components. The interior of the star cools through the emission of neutrinos, which is very sensitive 
to composition and possible superfluidity. The crucial question is whether or not the neutrino emission 
processes are rapid or slow, terms having historical significance. The surface temperature of the neutron 
star, apart from an initial period of several to one hundred years, the thermal timescale of the crust, 
is intimately coupled to that of the core. After the crust reaches thermal equilibrium with the core, a 
rapidly cooling core would be revealed through relatively low surface temperatures for a given age. It 
is now realized that the behavior of the symmetry energy of dense matter is crucial in determining the 
relative rate of neutrino cooling. This is because the symmetry energy not only controls the rate of cooling 
for a mixture of nucleons, but also determines at what densities exotic material, such as hyperons, Bose 
condensates, or deconfined quark matter appear, all of which could allow relatively rapid cooling even if 
neutrino processes involving nucleons were slow. 

There are other, less traditional, ways in which information about neutron stars can be gleaned. 
These involve neutrinos and gravity waves. Indeed, the fewer than two dozen neutrinos observed from SN 
1987A are only the harbinger of the thousands expected in current detectors, or millions from planned 
detectors, from a galactic supernova. These neutrinos are emitted in two phases. The first is a tremendous 
burst accompanying core bounce, with peak luminosities of nearly 1000 bethes s" 1 (1 bethe = 10 51 erg), 
exceeding the photon luminosity of all the stars in the visible universe during the same instant. The second 
is a continuous flux, quickly ramping down from the peak of the initial burst, from the deleptonizing 
proto-neutron star that will remain visible for perhaps fifty seconds. Nearly all (~ 99%) of the neutron 
star's binding energy of about 300 bethes is eventually radiated away in neutrinos of all flavors . The 
neutrino luminosities and the emission timescale are controlled by several factors including the total mass 
of the proto-neutron star and the opacity of neutrinos at supranuclear density, which depend on the star's 
composition and the EOS of strongly interacting dense matter. The details of the neutrino light curve and 
neutrino energy distributions could reveal details of their internal compositions. 

Last, but perhaps not least, prodiguous amounts of gravity waves are expected to be emitted in the 
death throes of binary mergers involving neutron stars. Compact binaries, containing neutron stars and/or 
black holes, continuously radiate gravitational waves which cause their orbits to decay. Observed neutron 
star binaries have times-to-merger ranging from about 85 million years to tens of billions of years. Current 
estimates imply that there are about 1 merger per 10,000 years in our galaxy, and detectors such as 
the advanced LIGO might observe 40-400 events per year [7|. The observed patterns of gravitational 
waves could yield information about the masses of the components as well as estimates of the radii of 
inspiralling neutron stars. Furthermore, as we discuss below, there should be a pronounced difference in 
the observed wave pattern between mergers involving normal and self-bound strange stars. In contrast to 
other observational differences between normal and strange stars, which could have multiple sources, the 
differences observed in gravitational radiation would be relatively unambiguous and therefore, possibly, 
unique. 

Complementary to observational constraints on neutron star structure are laboratory measurements 
which can restrict the range of nuclear parameters such as the incompressibility, the symmetry energy and 
its density dependence, nucleonic effective masses and nuclear specific heats, the high-density pressure- 
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density relation, and hyperon-nucleon couplings. We briefly review laboratory work involving dense matter 
and its relation to neutron star structure and compositions (we provide references to more detailed sum- 
maries). This experimental information ranges from measurements of nuclear masses (including high-N 
and low-Z nuclei in Radioactive Isotope Accelerators), nuclear charge radii, neutron skin thicknesses, fre- 
quencies of giant monopole and dipole resonances, and collective flow, multifragmentation and isospin 
diffusion information from heavy ion collisions. 



A. Personal note (James Lattimer) 

I would like to take this opportunity to add some personal recollections of interactions with Hans Bethe. 
I was introduced to Hans in the summer of 1978 when I was an Urbana postdoc visiting Gerry Brown 
in Copenhagen. Gerry and Hans were interested in my work with Geoff Ravenhall on finite temperature 
nuclear phase coexistence jl] and Hans had developed his "low entropy" thesis that was so valuable to 
understanding gravitational collapse supernovae. During this time, I was invited to become a co-author 
of the paper that we, at least, have since referred to as BBAL (9|. I was greatly honored, of course. But 
when I was given a draft of the paper and saw it began with the phrase "Massive stars live for eons and 
eons I had second thoughts. I conveyed my misgivings about the astronomy to Gerry, but he assured 
me that Hans was very happy to accept criticism. I agreed to collaborate and began rewriting, much to 
Gerry's chagrin. After my return to Urbana, and Hans to Ithaca, Gerry would distribute my rewrites to 
Hans and Hans' to me. This iteration did not seem to rapidly converge, the paper developed into a 46-page 
journal article with seven appendices, and did not get submitted until the following Spring. Gerry still 
blames me for the delay, probably because a competing letter (with not dissimilar conclusions) I wrote at 



Urbana with Don Lamb, Chris Pethick and Geoff Ravenhall [10j had appeared shortly beforehand. 

Another memorable occasion occurred in 1981, a year or two after I came to Stony Brook (partly, at 
least, on the back of BBAL). I had frequently cooked breakfast and dinner for Hans during summer trips 
to Copenhagen over the previous few years; he especially loved meat and always liked the rare steaks I 
cooked for him in Denmark (at least, what passed for steaks in Denmark). Gerry decided I should host a 
dinner for Hans at which all the graduate students in our group would attend. I decided to cook a leg of 
lamb on the barbecue. Around this time, Hans had become interested in neutrino diffusion in supernovae, 
and was intrigued by comments from Don Lamb that the boundary in collapsing material between trapped 
and freely-streaming neutrinos was particularly important in supernovae. Hans referred to this theory as 
the joint of Lamb. So, of course, when I took my leg of lamb off the grill and put it on the table before 
Hans, I pronounced it "the joint of lamb" . Hans, who loved puns, had a great laugh over this. 



B. Personal note (Madappa Prakash) 

I distinctly recall my first meeting with Hans Bethe. One morning during my early years as a postdoc 
at Stony Brook, Gerry Brown came into my office and said "Come and tell Hans what you're doing now; 
we're in the common room", and left. Expecting to see Hans Hansson, a friendly particle physicist, I 
grabbed my coffee cup and went into the common room. Sitting there was Hans Bethe, eating raisins one 
by one. Gerry said, "Go on, write something on the board and entertain us." I recall feeling very ill. I 
said, "Could I bring my notes?" Gerry said, "Ok, if you really need them." I ran into my room, grabbed 
my notes and returned to the common room where the two of them were waiting. I started writing on 
the board feeling more and more ill, but fortunately I was interrupted. Hans Bethe said, "Why are you 
keeping so many terms in the expansion?" (The physics issue was the influence of effective masses on 
the pressure of supernova matter.) I somehow found the courage to say, "I also have an exact and short 
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expression." Bethe seemed to like that and said "If you have to keep more than one term in the expansion, 
you're expanding about the wrong point." Later, Gerry told me that Hans liked what I was doing, which 
was very gratifying. The work I was reporting to Bethe was later written up (with Tom Ainsworth, Jean- 
Paul Blaizot and Herman Wolter as co-authors) for Gerry Brown's 60th Birthday Conference Proceedings 



"Windsurfing the Fermi Sea" [11 



My most memorable meeting with Bethe took place in Seattle during the "Supernova Physics" pro- 
gram, organized by Gerry Brown during the early years of the U.S. Institute for Nuclear Theory and held 
in a building some distance away from the old Physics Department. I was in awe that Bethe had solved 
two coupled differential equations to three decimal point accuracy using his slide rule, while others were 
still cajoling their computers. One afternoon, Gerry asked me to bring Bethe to the brown-bag lunch 
in the Physics Department, where Stan Woosley was to speak. We started walking toward the Physics 
building. Bethe walked slowly, so we were getting late. To make matters worse, he stopped now and then 
to pop raisins into his mouth. We arrived at a street crossing light, which turned green as we approached 
it. I hurried Bethe as best as I could, but when we were half-way across the lights turned red again. To 
my utter horror, the perpendicular traffic was accelerating toward us! I stood in the middle of the road 
with my arms widespread to stop the racing cars. Fortunately, Bethe ambled across to safety, while I was 
busy stopping impatient drivers in the middle of the road! Imagine the thoughts in my head! (I would 
be forever known as the man who got Bethe killed!) When we arrived at the lunch meeting, I was really 
angry at Gerry, who just smiled and said "Why did you think I asked you to bring him here?" 

I wish there were more people like Hans Bethe. He was very kind to me, always made me feel good, 
and I learned a whole lot from him. I am very fortunate to have known him. 



II. SOME MAXIMAL ATTRIBUTES OF NEUTRON STARS 

A. The maximum mass 

General relativity introduces a profound change to Newtonian hydrostatic equilibrium: the existence 
of a maximum mass. In general relativity, hydrostatic equilibrium is expressed by 

dp{r) G [p(r) + e(r)][M(r) + 4:7ir 3 p(r)/c 2 } dM(r) 2 

~dT~~~^ r{r-2GM{r)/c 2 ) ' ~dr~ ~ * p{T > ' ^ 



where e = pc is the total mass-energy density. Rhoades & Ruffini [5| derived an upper limit of 3.2 M Q 
to neutron star masses under the assumptions that the EOS (1) is less stiff than that of non-interacting 
degenerate neutron matter up to a fiducial mass (energy) density (pf)ef, which they chose to be p/ = 
€f/c 2 = 4.6 x 10 14 g cm -3 , and (2) is limited by causality at higher densities. In other words, they assumed 
the EOS 

P(e) = Pneutron(t), e < e fi P(<0 = C ~ e / + V f , e > e /> ( 2 ) 



where pf is the density at ef. Hartle & Sabbadini [12j later showed this limit almost perfectly scales with 
e/ such that 



M max = 4.2 Je f /e s M . (3) 
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Scaling results for "minimum period" EOSs 



Eq. 
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2.825 


0.007247 


4.476 x 10 15 


0.02193 


13.55 x 10 15 


3.026 


0.200 





2.87 


0.00725 


4.478 x 10 15 


0.0219 


13.54 x 10 15 


3.026 


0.203 



TABLE I: The quantities R 1 ^^ , anc ^ PcerS, sph re f er to the radius, mass and central mass density of the 

maximum mass configuration. P m in is the Keplerian rotational limit of the maximum mass configuration. The 
first line of results corresponds to Eq. ([6]); the second line is for Eq. (J7]). 



B. Maximum compactness 



Concerning the maximum compactness, Lindblom [13J determined that the maximum redshift, for 
P/>3xl0 14 g cm -3 and assuming causality, is 

1< 0.863, (4) 



'I - 2GM/Rc 2 
which is equivalent to 

R > 2.83GM/c 2 . (5) 

This result is relatively insensitive to the value of pf and is consistent with the empirical limit established 
by Glendenning (l4j . 

Koranda et al. [HI have given a very elegant and rigorous causal limit that is independent of any 
other assumptions about the EOS. They emphasized the special nature of the "minimum period" EOS (so 
named for reasons made evident below) 

p(e) = 0, e < e c ; p(e) = e - e c , e > e c . (6) 

This EOS contains a single parameter, e c or p c , and as a result the general relativistic structure equations, 
Eq. (CQ), then contain only one free parameter and one dimensional parameter, e c , which corresponds 
to the surface energy density. All properties of spherical (i.e. non-rotating) stars with this EOS must 
therefore scale as M oc R and e oc M~ 2 . These considerations must apply, in particular, to the maximum 
mass configuration, whose radius, mass and central density are R^ l h x ,Ml^ h lx and e™^? h (or p™^ sp/t ), 
respectively. Furthermore, this being the maximally stiff EOS at high density, and the minimally stiff 
EOS at low density, it will support the largest mass with the smallest radius. For a given value of e c , the 
maximum mass configuration obviously is the most compact configuration. Integration of the dimensionless 
structure equations then yields the scaling relations shown in the first line of Table d for the maximum 
mass configurations. 

Of course, the EOS in Eq. ([6]) is unrealistic in that the density at the surface of the star is finite, with 
a value generally well above the nuclear equilibrium density. A less rigorous, but more realistic, causally 
limited EOS consists of using a reasonably understood model (here chosen to be the FPS EOS [l6j]) for 
p(e) up to a matching energy density e m where the pressure is p m : 

pie) = p FPS {e), e < e m ; p(e) = p F ps(^m), e m < e < e c ; p(e) = e - e c + p m , e > e c ■ (7) 

Between e m and e c , the pressure is fixed, and above e c , the EOS is causal. Choosing a matching baryon 
number density n m = 0.1 fm -3 , one has for the FPS EOS e m = 1.262 x 10~ 4 km -2 and p m = 5.69 x 10~ 7 
km -2 (units where c = G = 1). The maximum mass configurations retain, to a high accuracy, the scalings 
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FIG. 1: The maximum mass - central density relation predicted by causality coupled with the Tolman VII and 
Tolman IV analytic GR solutions are compared with structure integration results for a variety of EOSs. NR 
refers to non-relativistic potential EOSs, R refers to relativistic field-theoretical EOSs, and Exotica refers to EOSs 
with considerable softening at high density due to kaon condensation or strange quark matter deconfinement. A 
possible redshift measurement of z = 0.35 is also shown. Figure taken from. Ref. 17]. 



of the EOS Eq. ([6]), but the integration constants are slightly altered, as shown in the second line of Table 
HI These stars have slightly larger maximum masses and radii for a given value of e c . Also note that the 
revised compactness limit is practically the same as established by Lindblom [13 . 



C. Maximum central density 

The "minimum period" EOS predicts the smallest radius possible for a given mass, with R > 
2.87GM/c 2 . Dimensionally, p oc MR~ 3 , so a plausible conjecture would be that this EOS also predicts 
the largest possible value of the central density for a given mass. However, this is not the case. Lattimer 
& Prakash [17| investigated the maximum density question, with an eye toward using measured masses to 
set upper limits to the density of static matter in the universe. Their empirical results demonstrate that 
the predicted value of M 2 p cent for the "minimum period" EOS is close to, but less than, the actual limiting 
value. They found that an empirical limit can be found by combining the causality limit R > 2.87GM/c 2 
with the central density-mass relation implied by the Tolman VII analytic solution 18| of the general 
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relativistic structure equations. The Tolman VII solution has the explicit energy density-radius relation 

(8) 

which, combined with the causality limit, results in 

« 15 /M\ 3 15 / c 2 \ 3 „ 9 , 

This relation, and results for various EOSs, are displayed in Fig. [TJ From Tabled note that the "minimal 
period" EOS predicts a maximum central density with a coefficient of 13.5 x 10 15 in Eq. (Q, about 10% 
smaller than the empirical limit. 

It should be emphasized that Eq. (jHD represents the upper limit to the central density of the measured 
star. Since the maximum neutron star mass must be larger than that of any measured star, and the central 
density increases with mass for a given EOS, the central density of the actual maximum mass star must be 
smaller than the value given by Eq. (Q using the largest measured mass. With a mass measurement of 2.1 
M , for example, the limiting mass density is 3.4 x 10 15 g cm -3 . If the maximum mass was in fact about 
10% larger, the limiting mass density becomes 2.8 x 10 15 g cm -3 . Note that this mass density corresponds 
to a baryon density of only about 7n s . This could be small enough to call into question the applicability 
of perturbative QCD in invoking the presence of a quark phase in neutron star cores. It must be stressed, 
however, that precocious perturbativeness may indeed occur in the quark matter sector. 




D. Maximum spin rate 

The compactness limits above are also intimately connected to the Keplerian, or mass-shedding, rota- 
tional limit obtained when the equatorial surface velocity equals the orbital speed just above the surface. 
The most compact star, assuming uniform rotation, has the highest rotational frequency. For a uniform 
rigid sphere of mass M and radius R, the mass-shedding limit in Newtonian gravity is 



It is interesting and fortunate that the minimum spin period in fully relativistic calculations 
employing realistic hadronic EOSs is given to a good approximation by a similar formula 

P min ~ 0.83 — — 7 — ms. (11) 

\ M ™h ) V 10km / 

This spin period obtains for the maximum mass configuration, which has the smallest radius for hadronic 
EOSs. This formula takes into account not only general relativity, but also the deformation of rotating 
stars, yet is still expressed in terms of the non-rotating values of the maximum mass and radius. Therefore, 



the same scaling P m in oc JR 3 /M exists for both Newtonian and GR gravitation. 



The same scaling exists for the "minimum period" EOS [151 ] : the scaling R 1 ^ oc M™% x for maximum 
mass stars implies that Pmin oc M^ x . Koranda et al. [Hj determined the proportionality constant 0.200 
ms M" 1 for this case (see Table [I]). This translates into a coefficient of approximately 0.74 in Eq. pip . 
The ultra-compactness of stars constructed with the "minimum period" EOS accounts for the considerably 
smaller coefficient compared to that originally calculated for more realistic EOSs. 
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FIG. 2: Mass-radius trajectories for typical EOSs (see [6J for notation) are shown as black curves. Green 
curves (SQM1, SQM3) are self-bound quark stars. Orange lines are contours of radiation radius, i?oo = 
R/\/l - 2GM/Rc 2 . The dark blue region is excluded by the GR constraint R > 2GM/c 2 , the light blue re- 
gion is excluded by the finite pressure constraint R > (9/4)GM/c 2 , and the green region is excluded by causality, 
R > 2.9GM/c 2 . The light green region shows the region R > R ma x excluded by the 716 Hz pulsar J1748-2446ad 



22| using Eq. (|12J) . The upper red dashed curve is the corresponding rotational limit for the 1122 Hz X-ray source 
XTE J1739-285 [23(| ; the lower blue dashed curve is the rogorous causal limit using the coefficient 0.74 ms in Eq. 
121). 



The relation in Eq. (lllj) applies only to the maximum mass configuration. Lattimer & Prakash [2 11 ] 
empirically found a somewhat more useful result which applies to an arbitrary neutron star mass, so long 
as that mass is not close to the maximum mass: 

P m ^ (0.96 ± 0.03) (J£) (J?£) ms . (12) 

In this equation, M sp h and R sp h refer to the non-rotating mass and radius. Inasmuch as an observed 
neutron star is likely to be at least slightly smaller than the maximum mass, Eq. f|T2l) can be used to limit 
masses and radii for observed stars: 

/1000 Hz \ 2/3 fM svh \ 1/3 , , s 

R, ph< l 0A (-^) km, (13) 



12 



where the spin frequency v — 1/P. 

The most rapidly rotating pulsar is PSR J1748-2446ad with a spin rate of 716 Hz [22J. With this 
value, Eq. (fT3|) suggests that for a 1.4 M star, the non- rotating radius would be limited by R < 14.3 km 
as shown in Fig. [2J However, there is some liklihood that this star spun up through the accretion of a 
few tenths of a solar mass. If this star is, for example, 1.7 M Q , the upper limit to the non-rotating radius 
would be 15.3 km. Unfortunately, neither limit is very restrictive at present. 

Recently, however, an 1122 Hz X-ray burst oscillation from the neutron star X-ray transient XTE 
J1739-285 has been reported [23]. The stability of this oscillation frequency strongly suggests that it is 
the spin rate of the neutron star. If true, this sets relatively stringent limits to the radius that this star 
would have had if it was not rotating. According to the empirical relation Eq. ffl~3]) . the maximum radii 
are 10.8 and 12.2 km, assuming the mass of this star to be 1.4 and 2.0 M , respectively. Generally, these 
limits apply for stars not in the immediate vicinity of their maximum mass or for quark matter stars. As 
shown in Fig. [2] this constraint would rule out a number of EOSs. The rigorous causal limit, from Table 
[U yields 12.6 and 14.2 km, respectively. It will obviously be important to confirm this observation. 



III. RECENT MASS MEASUREMENTS AND THEIR IMPLICATIONS 

Several recent observations of neutron stars have direct bearing on the determination of the maximum 
mass. The most accurately measured masses are from timing observations of the radio binary pulsars. 
As shown in Fig. [31 which is compilation of the measured neutron star masses as of November 2006, 
observations include pulsars orbiting another neutron star, a white dwarf or a main-sequence star. The 
compact nature of several binary pulsars permits detection of relativistic effects, such as Shapiro delay 
or orbit shrinkage due to gravitational radiation reaction, which constrains the inclination angle and 
allows the measurement of each mass in the binary. A sufficiently well-observed system can have masses 
determined to impressive accuracy. The textbook case is the binary pulsar PSR 1913+16, in which the 
masses are 1.3867 ± 0.0002 and 1.4414 ± 0.0002 M G , respectively 0. 

One significant development concerns mass determinations in binaries with white dwarf companions, 
which show a broader range of neutron star masses than binary neutron star pulsars. Perhaps a rather 
narrow set of evolutionary circumstances conspire to form double neutron star binaries, leading to a 



restricted range of neutron star masses [53j. This restriction is likely relaxed for other neutron star 
binaries. Evidence is accumulating that a few of the white dwarf binaries may contain neutron stars 
larger than the canonical 1.4 M value, including the intriguing case [45[ of PSR J0751+1807 in which the 
estimated mass with la error bars is 2.1 ± 0.2 M . In addition, to 95% confidence, one of the two pulsars 



Ter 5 I and J has a reported mass larger than 1.68 M (43 . 

Whereas the observed simple mean mass of neutron stars with white dwarf companions exceeds those 
with neutron star companions by 0.25 M , the weighted means of the two groups are virtually the same. 
The 2.1 M neutron star, PSR J0751 + 1807, is about 4a from the canonical value of 1.4 M . It is 
furthermore the case that the la errors of all but two systems extend into the range below 1.45 M Q , 
so caution should be exercised before concluding that firm evidence of large neutron star masses exists. 
Continued observations, which will reduce the observational errors, are necessary to clarify this situation. 

Masses can also be estimated for another handful of binaries which contain an accreting neutron star 
emitting x-rays, as shown in Fig. [31 Some of these systems are characterized by relatively large masses, 
but the estimated errors are also large. The system of Vela X-l is noteworthy because its lower mass limit 
(1.6 to 1.7M ) is at least mildly constrained by geometry 26l . 



Raising the limit for the neutron star maximum mass could eliminate entire families of EOSs, especially 
those in which substantial softening begins around 2 to 3n s . This could be extremely significant, since 
exotica (hyperons, Bose condensates, or quarks) generally reduce the maximum mass appreciably. 
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FIG. 3: Measured and estimated masses of neutron stars in radio binary pulsars (gold, silver and blue regions) and 
in x-ray accreting binaries (green). For each region, simple averages are shown as dotted lines; weighted averages 



are shown as dashed lines. The labels (a) = 2jj through (C) = 52( are references cited in the bibliography. For 



the stars with references z-C, a lower limit to the pulsar mass of 1 M Q was assumed. 



Assuming that the hyperon-nucleon couplings are comparable to the nucleon-nucleon cou plings typi 



55. 56. 57 



cally results in the appearance of A and E~ hyperons around 2 to 3 n s in neutron star matter [5 
In beta equilibrated neutron star matter, the various chemical chemical potentals satisfy the relations 
— H P = He = A*s- and fi n = fi\. As a consequence, the proton fraction in such matter is quite small, 
of order 5-10%. Little is known about the symmetry dependence of the hyperon-nucleon couplings as 
these couplings are chiefly determined from hyperon binding energies in more or less symmetric nuclei. If 
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hyperons indeed appear at as low a density as 2-3 n s , the maximum neutron star mass becomes relatively 
small, typically less than 1.6 M Q jH^j. 

The suggestion of Kaplan and Nelson [58[ that, above some critical density, the preferred state of matter 
mig ht contain a Bose-Einstein condensate of negatively charged kaons has been examined extensively 



62, |63j, |64j, |65|, |66|, |67|, |68j. The astrophysical consequences have been explored in some detail in Ref. 



68] . The physics is that in medium, the strong attraction between K~ mesons and baryons increases with 
density and lowers the energy of the zero-momentum state. A condensate forms when this this energy 
becomes equal to the kaon chemical potential, /i which is related to the electron and nucleon chemical 
potentials by /i = p n — fi p = fi e = p^ due to chemical equilibrium in the various reactions. Typically, the 
critical density for condensation (which depends primarily on the symmetry energy of nucleonic matter) is 
~ (3 — 4)n s , although it is model and parameter dependent. Relative to matter without a kaon-condensed 
state and depending upon the models employed, maximum masses only as high as 1.5 to 1.6 M Q can be 
obtained with kaon condensation. 

If a different form of strangeness can appear prior to hyperons or Bose condensates, for example, 
deconfmed u, d, and s quark matter, maximum masses up to approximately 2.0 M Q are possible (69[, but 
only with special fine-tuning of the nucleonic and quark matter parameters. Therefore a confirmation of 
a neutron star mass in excess of 2 M would be especially interesting. 

A parameter of increasing interest, as the database of neutron star masses expands, concerns the min- 
imum mass of neutron stars. It is believed that neutron stars are created in the aftermath of gravitational 
collapse supernovae and must therefore pass through a proto-neutron star state in which a large number 
of neutrinos are trapped and a moderate entropy exists [1]. For such a configuration, the range of stable 
masses is much narrower than for cold, catalyzed systems. The minimum mass of a cold star is about 0.09 
M Q , only slightly sensitive to the EOS around p s . A star with a trapped lepton fraction Yl ~ 0.35 — 0.4 
and an entropy per baryon s ~ 1 — 2 has a minimum mass of about 0.85 M [70j] for stability. The smallest 
reliably estimated neutron star mass is the companion of the binary pulsar J1756-2251 whose mass is 
1.18 ± 0.02 M Q. It will be interesting to see if smaller neutron stars are found. 



IV. THE RADIUS CONSTRAINT 

As previously mentioned, the radius of a neutron star is primarily determined by the density depen- 
dence of the symmetry energy. This connection arises through the relation between the neutron star radius 
and the internal pressure of matter at intermediate densities (1.5n s < n < 2 — 3 n s ). That such a relation 
exists can be seen most clearly by considering Newtonian polytropes. For the pressure-density relation 

p = Kp 1+1 l\ (14) 

where K is a constant and n is the polytropic index, hydrostatic equilibrium implies that 

R oc K n/{3 ~ n) M {1 - n)/{3 - n) . (15) 

Realistic EOSs typically have n ~ 1, but K is uncertain by a factor of 5 or 6. For some average density 
p*, in the vicinity of (1 — 2)mbn s , suppose the pressure is p*. For the case n — 1, one therefore has 

Rtxpl^p^M . (16) 

The independence of R from M is a characteristic result of the mass-radius trajectories for EOSs without 
extreme low-density softening, as seen in Fig. [2j Equation ffTBI) suggests that the radius scales with the 
square root of the fiducial pressure p*. However, general relativity plays an important role, and with its 
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FIG. 4: Empirical demonstration of the constancy Rp* , for 1 M Q (upper panel) and 1.4 Mq (lowerpanel) 
stars. For each mass, 3 fiducial number densities are selected. Figure and EOS labels are from reference [6(. 



inclusion, the scaling is such that R oc pj [6|. This is an empirical result, and is graphically illustrated 
in Fig. HI It is important to note that it applies to nearly all neutron star EOSs, not just those that have 
a portion of the M — R curve with vertical slope. The exceptions to the R oc pV 4 rule are the EOSs for 
strange quark stars or EOSs with extreme softening near p s . The correlation is more robust if the fiducial 
density is taken to be 2n s rather than n s . 

Using the only known analytic solution of Einstein's equations with vanishing surface energy density 
and an explicit pressure - density relation, this scaling can be analytically demonstrated. This solution, 
due to Buchdahl [7lj], assumes the EOS 

e = 12y/p~p - 5p, (17) 

where p* is a parameter. For low densities, when p « p*, one sees that this EOS is that of an n = 1 
polytrope, so this solution is a reasonable approximation for a neutron star. For this solution, the radius 
is an implicit function of M and p*: 

^C-^feFl' (18) 
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where ft = GM/Rc 2 . In turn, these relations lead to 

dlnR 



dlnp 



n,M 



l (l-/3)(l-2/?)6-5 v ^M 



Choosing for typical values M = 1.4 M Q ,R = 13 km, and e = 1.5m;,n s = 2.97 x 10 4 km 2 , one finds 
ft = 0.318 and, using Eq. dTTJ), p/p* = 0.141. Substitution in Eq. (JUD yields GHniZ/dlnp = 0.230, quite 
close to the empirical result of about 1/4. 

This correlation is significant because the pressure of degenerate neutron-star matter near the nuclear 
saturation density n s is, in large part, determined by the symmetry properties of the EOS. For the present 
discussion, we introduce the incompressibility K and the skewness K', and expand the nucleonic energy 
per particle about its values at n s and x = 1/2, where x is the proton fraction: 

£(n,z) = -16 + -(l--) -2,)'.... (20) 

Here, E sym is the symmetry energy function, approximately the energy difference at a given density between 
symmetric and pure neutron matter. The symmetry energy parameter is defined as S v = E sym (n s ). For 
the total energy, leptonic contributions (mainly from electrons as that from muons is small in the vicinity 
of the nuclear equilibrium density) 

E e = (3/4)hcx(3ir 2 nx A ) 1/3 (21) 

must be added to E(n,x). Because catalyzed matter in neutron stars is in beta equilibrium, i.e., /i e = 
yU n — yLtp = —dE/dx, the equilibrium proton fraction at n s is 

x s ~ (?,'K 2 n s )- l (AS v /hcY ~ 0.04 . (22) 

This small value of x s enables the pressure at n s to be expressed as 

%x s )=( / 1% y cKl)1 | = n s (l-2x s )[n s S' v (l-2x s ) + S v x s ]^n 2 s S' v , (23) 

V / n s ,x s 

where S' v = (dE sym /dn) ng . The pressure depends primarily upon S' v , because the terms proportional to 
x s are relatively small. The equilibrium pressure at moderately larger densities is similarly insensitive to 
K and K' . Experimental constraints to the compression modulus K, most importantly from analyses of 



giant monopole resonances give K = 220 MeV [72J. The skewness parameter K' has been estimated to lie 



in the range (1780-2380) MeV |73j. Evaluating the pressure for n = 1.5n s 



P(1.5n s ) = 2.2571^/18 - /i'/216 + n s (l - 2x) 2 (dE sym /dn) 1 . 5ris ] . (24) 

Note that the contributions from K and K' largely cancel leading to the result that the pressure is 
dominated by the term involving the density derivative of the symmetry energy. 

The density dependence of the symmetry energy is constrained by laboratory data such as nuclear 
mass fits, giant dipole resonances, and neutron skin thicknesses of neutron-rich nuclei. However, significant 
difficulties are encountered in obtaining the quantity S' v in the necessary supra-nuclear density range. First 
of all, nuclei sample nucleonic matter in the density range n < n s . However, to determine neutron star radii, 
we need the symmetry energy in the range n s to 2n s . Secondly, nuclear mass fits and giant dipole resonance 
data predict a strong correlation between the volume symmetry energy S v and the surface symmetry energy 
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S s parameters in a liquid-drop sense. Neutron skin thickness is predicted to be proportional to the ratio 
S s / S v . If the nuclear energy per baryon can be expressed as in Eq. (1201) . a simplified Thomas- Fermi model 
for the nucleus shows that S s / S v and the neutron skin thickness are proportional to the integral quantity 



Refs. U4, 75 



(S v /E, 



sym 



n 



zdn . 



y/n(E(n, 1/2) + 16) 

Although S' v is constrained by this relation, it cannot be uniquely determined from it. 



(25) 



V. NEUTRON STAR CRUSTS AND THEIR CONSTRAINTS 

There are a number of observations that pertain to properties of the neutron star crust, which specifi- 
cally is the region between the nuclear surface and the phase transition separating uniform nuclear matter 
and matter containing nuclei. The phase transition occurs at about 1/3 to 1/2 p s , or about 10 14 g cm -3 , 
depending on the nuclear compression modulus and the density dependence of the nuclear symmetry en- 
ergy. For example, if pulsar glitches are due to interactions between the neutron superfluid confined to 
the crust and the bulk of the star, the rate of angular momentum transfer can be related to the fraction 
of the moment of inertia of the star which resides in the crust [76]. Quasi-periodic oscillations observed 
in X-ray emission following X-ray bursts on neutron stars are likely related to fundamental and overtone 



vibrational frequencies of neutron star crusts [77] . Finally, the cooling observed over the first several years 
following superbursts from neutron stars or giant flares from magnetars is likely due to the cooling of the 
crust and can be a sensitive indicator of its extent [781 . 



A. Theoretical considerations 

It is possible to relate both the relative thickness of the crust, A/R, and the relative fraction of the 
moment of inertia contained in the crust A///, to the mass, radius and a single parameter of the core-crust 
interface which depends on the EOS. To see this, first consider the crust thickness A = R — R t where R t is 
the radius of the core-crust interface. Begin with hydrostatic equilibrium in the crust, setting M(r) = M, 
ignoring P relative to both e and to Mc 2 /47i\R 3 , and ignoring the internal energy per baryon compared to 
rribC 2 . Then, since p ~ m^n, 

Mr) = dp = GM ^ ^ 

m^n nib r 2 — 2GMr/c 2 ' 

where p — p n is the baryon chemical potential in beta equilibrium. Integrating from the base of the crust 
(values of physical variables here are denoted by the subscipt t) to the surface (denoted by subscript 0), 
Eq. (1261) becomes 

p t -po _ l^ r t (R-2GM/c 2 ) 
m b c 2 2 n R(r t — 2GM/c 2 ) ' 1 ' 

Note that r§ = R and r t = R — A. Approximately, fiQ ~ — 9 MeV. Upon defining 

n _ e 2( Mt - M0 )/m 6 c 2 ) ( 2g ) 



Eq. (}2~7l) can be manipulated to yield 
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In general, H-l « 1, whence Eq. ([29]) implies that A oc (1 - 2/3)R 2 /M. 

In a similar way, one can determine the crustal moment of inertia fraction. One begins with the 
definition of the moment of inertia in general relativity [79| : 



2r 2 r R di(r) Kir r R 

I r 3 u ^^l dr = lL r\p{r)+p{r)/c 2 )e x{r) j{r)u{r)dr, (30) 



3G Jo dr 3c 2 Jo 

where j(r) = e ~>M+ A M)/ 2 and oj{r) is the rotational drag. The metric functions u(r) and A(r) satisfy 



<W) = 2C M(r)+p(r)/c 2 ^ e _ A(r) = x _ 2GM{r) 
dr r(r — 2GM(r)/c 2 ) ' rc 2 



(31) 



whereas the rotational drag satisfies 



^ m ^y. irMr) m. (32) 

The relevant boundary conditions are 

■mo) = dm = *<o) = A4o) = o 

<ir c/r c?r 

Substitution of Eqs. ([32]) and ([33]) into Eq. ([30]) yields 

(34, 

While these equations are straightforward to solve for a given EOS, it is useful to employ a fit estab- 
lished by Lattimer & Schutz 8(| which is valid for realistic hadronic EOSs that permit maximum masses 
greater than about 1.6 M Q : 

/ ~ (0.237 ± 0.008) MR 2 (1 + 2.84/3 + 18.9/? 4 ) M km 2 = aMR 2 f((3) , (35) 

which defines a and f(/3) for future reference. 

One can estimate the moment of inertia connected with the crust by employing the right-hand form 
of Eq. ([30]) together with hydrostatic equilibrium, Eq. (JTJ, and the approximations 4Trr 3 p(r) << 
M(r)c 2 , M(r) ~ M,j(r)u(r) ~ uj(R), to find 

A/= /^ A r 4 (p + p/c 2 )e- A Mj(r)^(r)cir ~ /f A r >- ( 36 ) 

If one assumes that p oc p 7 in the crust, and that p(i? — A) = p t , the integral / r 6 dp is approximately 



rR ( Qry 

/ r 6 dp ~ /2 6 pi exp 

Jfl-A \ 27 — 1 



67 A s 
R 



(37) 



Combining Eqs. ( [38]) . ( [33]) and ( [35]) . one finds 

AI 8tt i^pt / 67 A 



3c2M2L a- 1 /(«- 1 -2/3]exp . (« 
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Considering the importance of the location of the core-crust interface, it is worthwhile analyzing how 
this depends upon EOS assumptions. A heuristic determination of H, can be made as follows. The 
core-crust interface corresponds to the phase transition between nuclei and uniform matter. The uniform 
matter is nearly pure neutron matter, with a proton fraction of just a few percent determined by the 
condition of beta equilibrium. Ignoring the finite size effects due to surface and Coulomb energies of 
nuclei, one can consider the critical density for which uniform npe matter becomes unstable to separation 
into two coexisting phases (one corresponding to nuclei, the other to a nucleonic sea) [8l| . This is a first 
approximation to the actual case, in which surface and Coulomb effects are not negligible and in addition 
the dense phase disappears when its volume fraction is finite (82|. The weak interactions conserve both 
baryon number and charge. At zero temperature, the first law of thermodynamics can be written as 



du = —pdv — fidq . 



(39) 



where u is the internal energy per baryon, v = 1/n, q = x — Y e is the charge density, x is the proton 
fraction in the baryons, Y e is the electron concentration, and fi = a n — fi p = a e from beta equilibrium. 
Stability of the uniform phase requires that e(v, q) is convex 83j, or 

( dp\ ( dp\ I dq 



dv 



dq 



V \ / [L 



> o, 




(40) 



We assume baryonic and electronic energies of the form Eqs. (1201) and (I2TI) . which we rewrite as 



E N ~ V{n) + E sym {n){l - 2xf , E e = -Y e fi . (41) 

The beta equilibrium condition is 

fi = AE sym (n){\ - 2x) = hc{3ir 2 nx) 1/3 , (42) 
where charge neutrality x = Y e is used. The conditions in Eq. (1401) therefore become equivalent to 



2 d 2 V dV , N2 
n 2 — + 2n— + (l-2x 2 
dn z dn 



n 



2 d E S y m dE S y m j-i— i i dE S y m 

+ ln — A ^sym n — A 

dn y \ dn 



dn 2 



> 



'dq\ 1 



31; 
fi 



> 



(43) 



The second of these is always true, so the first one determines stability. 

As an example, consider the following simple density dependences for V and E sym : 



K f n 
V ^-- lQ + ~A l -n, 



E sym {n) = T k ( -) 2/3 + (S v - T k ) (- 



n. 



n. 



(44) 



where the kinetic contribution to the symmetry energy is T k — 13 MeV. A convenient parameter describing 
the density dependence is S' v , defined by 



n f 



iS,, 



'dlnE 



sym 



dlnn 



2T fc /3 + z(^-T fc ) 



(45) 
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FIG. 5: Contours of enthalpy Tit, transition density nt/n s , and pressure pt are displayed as functions of S v and 
the logarithmic derivative of E sym at n s for the simple model described by Eq. (j44|) . 



Applying the stability condition Eq. (143]) determines the phase boundary density n t . For reasonable ranges 
of S v and S' v , one finds 0.5 < n t /n s < 0.7 as seen in Fig. [51 The quantity H is related to the neutron 
chemical potential, 

finint) = -16 + ^ (l - ^) (l - 3^) + E sym {n t ){l - 4s?) + (n ^ (n) ) (1 - 2x t f , (46) 



by Eq. ( |28|) . The pressure at the boundary, is given by 



+ n , (1 _ 2ll 

9 n, Vn 9 



(47) 



The quantities 7i and p t at the phase boundary are also shown in Fig. [5] for this simple model, using 
n s = 0.16 fm -3 and K = 225 MeV. The boundary density n t and pt are much more sensitive to the density 
dependence of the symmetry energy than to S v or K. The general traits found in this simple model apply 
to realistic EOSs. In particular, for realistic EOSs 0], 

1.04 <H< 1.07; 0.20 MeV fm" 3 < p t < 0.65 MeV fm" 3 . (48) 

These ranges are measures of the current uncertainty in the density dependence of the symmetry energy. 
The range of values for Tit found for the simple model (and shown in Fig. [5]) is very similar. However, 
the realistic values for p t are somewhat smaller than predicted by the simple model, due to the neglect of 
finite-size effects. Finite-size effects are much less significant for the neutron chemical potential. 
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FIG. 6: The same as Fig. ([2]), but with additional contours of fixed redshift z = 0.35 (dot-dashed blue line), 
crustal fraction of moment of inertia AI/I = 0.14 (dashed blue line), and QPO constraints (red solid and dashed 
lines). For the Vela pulsar, the AI/I constraint and the assumption that the pressure at the core-crust interface 



satisfies pt < 0.65 MeV fm 3 implies that allowed masses and radii lie to the right of the dashed line 761 ] . 

B. Observational application 1: Pulsar glitches 

Pulsar glitches, which are sudden discontinuties in the spin-down of pulsars, seem to involve the transfer 
of angular momentum from an isolated component to the entire star. A leading model for glitches supposes 
that the isolated component consists of superfmid neutrons in the crust j84j. As the star spins down due 
to pulsar electromagnetic emissions, the superfmid component's angular velocity differs to a greater and 
greater extent from that of the non-superfmid component. When this differential becomes large enough, 
angular momentum is suddenly transferred from the superfluid component to the non-superfmid part 
through vortex unpinning. In the case of the Vela pulsar, some 30 years of observations indicate a steady 
overall angular momentum transfer rate that indicates at least 1.4 of the total moment of inertia of the 
star is involved with the isolated component 76[. Eq. (1361) indicats that the crustal fraction of the moment 



of inertia is most sensitive to the pressure at the transition from the crust to the core. The location of the 
transition and the pressure there are sensitive to, among other nuclear parameters, the density dependence 
of the nuclear symmetry energy. Using the upper limit of pt from Eq. PHI) in Eq. fl38|) allows one to set a 
minimum R for a given M for Vela U& which is displayed in Fig. [6] and is approximately: 



M 

i?>3.6 + 3.9 km. (49) 

" M 

Better understanding of the nuclear symmetry energy would permit a lower upper limit of p t to be deter- 
mined, which could tighten this constraint. 
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FIG. 7: Mass-radius diagram showing the constraints from neutron star seismology from the soft gamma-ray 
repeater SGR 1806-20. 7i values in the range 1.04 (lower boundaries of permitted regions) to 1.07 (upper bound- 
aries) were assumed. A few representative M — R curves from Fig. [2] are shown for reference. The straight lines 
follow the values of (3(7i = 1.04) = 0.116 and [3(71 = 1.07) = 0.165, determined from the observed frequencies, 
and show how the (M, R) values change as the shear velocities are varied. 

It should be pointed out that the standard picture of the outer core, in which superfluid neutrons co- 
exist with type II superconducting protons, has been challenged by Link jiBj. He notes that observations 
of the long-term precession in several pulsars indicates an incompatibility with creeping neutron vortices, 
which predict too large precession frequencies. Furthermore, the observations seem to require that the 
vortices of the inner crust must be able to move with little dissipation with respect to the solid, which 
could be incompatible with a neutron superfluid in the crust in precessing neutron stars. Nevertheless, 
since there is no evidence as yet that the Vela pulsar is precessing, Link's argument may not apply to this 
pulsar. Why some stars show precession and others do not is an unsolved problem, and until this is better 
understood, radius limits based upon glitches remain tentative. 



C. Observational application 2: Neutron star oscillations 



Another kind of limitation on radii can be obtained from observations of quasi-periodic oscillations in 



the X-ray emissions followi ng g iant flares in three soft gamma-ray repeaters, SGR 0526-66 
20 [8711 and SGT 1900+14 [88(1 , presumably due to torsional vibrations of the star's crust 



86] 
71 



SGR 1806- 
Neglecting 



and SGT 1900+14 [ 

magnetic fields and the effects of superfluids in the crust, the frequencies of the fundamental mode and 
overtones have been estimated [89j to be 



UJ 



n=0,£=2 



2v t 



RRt 



u n >o 



(50) 
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where v and A refer to the average values of the metric functions, and v r and v t are the average radial and 
transverse shear speeds, respectively. Higher order modes (£ > 2 in the case n = and I > 1 in the case 
n > 0) are ignored for a simple analysis. To a good approximation, one may assume e u = e~ A = 1 — 2/3. 
In the isotropic limit, v r = v t and we assume this; the shear velocities are of order 10 8 cm s _1 . Since the 
crust thickness A is small compared to R one has, roughly, tu n =o -R^, 1 , where Roo = R/y/l- — 2(3 is the 
so-called radiation radius (see Wip . Also, u> n>0 oc (1 — 2(3)/ A oc R 2 /M. 



Following the general approach of reference [89J, in which this simple analysis is compared to more 



sophisticated numerical models, one can write, using Eq. (l2"9j) for A/R, the frequencies (/ = lj/2tt) 



/km\ H - 1 + 213 , , 

f n>0 ~ 1170n f — J H _ l P - Hz . (51) 

If more than one frequency can be measured, these relations can be used to uniquely identify R and (3, 
that is, R and M, modulo uncertainties in Ti. 

The chief model dependence of the above approach lies in the velocities v r and v t . In the above, they 
are taken to be equal and are fitted to numerical simulations ignoring magnetic fields and superfluidity. 
In a more realistic analysis, they could differ from these values. However, if v r ~ Vt, Eqs. (I5ip indicates 
that (3(7i) can be found by eliminating R: it is the solution of a quadratic equation. 

The soft gamma-ray repeater SGR 1806-20 has quasi-periodic oscillations which correspond to the 
n = 0,£ = 2,6,10 and n = 1,£ = 1 modes [89| . The identification of the larger I modes is more 
speculative, but the frequencies f n =o,t=2 — 29 Hz and f n =i,£=i — 626.5 Hz are claimed. Fig. [7| shows the 
allowed regions in the mass-radius plane for the two modes, and their overlap, assuming 1.04 < H, < 1.07. 
It is clear that the overtones are much more sensitive to assumptions about Ti, (an EOS attribute) than 
is the fundamental mode. The two values of (3 associated with these values of 7i are (3{7i = 1.04) = 0.116 
and f3(H = 1.07) = 0.165, respectively. Because lines of constant (3 are approximately parallel to the 
overtone constraint, but perpendicular to the fundamental mode constraint, uncertainties in the velocity 
would be reflected as relatively larger uncertainties in the latter than the former. 



D. Observational application 3: Thermal relaxation of the crust 

After the formation of a proto-neutron star in the aftermath of a gravitational collapse supernova 
event, cooling of the core proceeds via neutrino emission, leading to temperatures of order 10 9 K within 
years (see WHip . During this period, the star is not in complete thermal equilibrium due to the relatively 
long thermal relaxation time of the crust, which is expected to be of order 10 to 100 years [9(| • The 
excess heat of the crust is preferentially transported by electron conduction into the star's isothermal core, 
dominating radiation by thermal emission from the surface. However, after a time which we will denote 
by t w , the surface will finally come into thermal equilibrium with the core, causing a reduction in the 
surface temperature. Depending upon whether or not the core neutrino emission processes are rapid or 
conventional (see the terminology in §VI Bp . the drop in surface temperature will be a factor of 2 (and 
relatively slow) to 10 (and relatively abrupt). Follo wing this reduction, the surface temperature stabilizes 
at a value about 100 times the core temperature [9l|. The time tn is simply related to the thermal 
conductivity and specific heat of the crust, its thickness, and redshift 9(J 92], so that its observation could 
be used to constrain neutron star structure and composition. 

In addition, some accreting neutron stars are observed to undergo long-duration X-ray bursts, some- 
times called superbursts, which appear to be powered by thermonuclear He fusion. During the burst, some 
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of the heat is transported into the interior, warming the crust. Observations of quiescent thermal emission 
following the bursts would therefore probe the thermal relaxation timescale of the crust as well as the 



rate of core neutrino emission [78j. It is therefore useful to investigate the relation between the thermal 
relaxation time and the structure of neutron stars. 

To begin, it is necessary to supplement the hydrostatic equilibrium equations, Eq. (TjQ), with those of 
relativistic heat transport, 

,dTe u ' 2 



—Anr Ke 



2 ^-( A -^)/ 2 ' 



dr ' 

%T = -4-V^C,f , (52) 

where L u (r) is the local neutrino luminosity, Cy is the specific heat capacity at fixed volume, and K is the 
thermal conductivity. Explicit neutrino emission is neglected in Eq. ( 1521) . Following thermal relaxation, 
the redshifted temperature T(r)e u ^ r ^ 2 becomes constant within the star. Applying these equations to the 

crust, and assuming A/R « 1 so that e u ~ e~ A ~ \A — 2GM/Rc 2 throughout the crust, and defining 
x — (R — r) / A, one finds 



r.\ 2 a 



OX 



Cve 3X/2 ^. (53) 



.Ay dx 

If it is assumed that the conductivity K and specific heat Cy can each be written as separable functions of 
density and temperature, then one can make the same assumption for T: T(r,t) = T ip(r)<p(t). Generally, 
Cy oc T, and K oc T° or T depending on density. The temporal dependence in these two cases is then 



either linear or quadratic: = 1 — t/t w or = Jl — t/t w , where the thermal lifetime of the crust is 

t w oc A 2 e 3A/2 , (54) 

with the proportionality constant scaling as Cy/K. This behavior is borne out in detailed numerical 
simulations of neutron star cooling 9(| in which, for superfluid crusts, 



t, n ~ 5 



f A\ 2 f 2GM\~ 3/2 
T- l--5^r 3^ (55) 



\km) V Rc 2 J 



For non-superfluid crusts, the cooling times are about a factor of 3.5 times larger. 

Although no newly-formed neutron stars have ever been observed so young as to show this surface 
temperature drop, it has been suggested that the resaon magnetars (highly magnetized neutron stars) have 
abnormally large surface temperatures is that they have heat sources (presumably powered by magnetic 
field decay) located in a thin layer in the outer crust [93|]. If the heat source was located deeper, the energy 
would be radiated efficiently by neutrinos and could not heat the crust. This heat source warms up the 
crust and produces an inhomogeneous temperature distribution within the star, with parts of the crust 
much warmer than the deep crust and core. The thermalization times calculated above are consistent with 
observations, but tell us more about the heat source than about structural constraints. 

This theory can be applied to superbursts or other cases in which the crust is heated for a long-enough 
duration that the crust thermalizes throughout. For example, observations of the soft X-ray transient 1H 
1905+000, which underwent a burst about 20 years ago, indicate that the crust cooling times are somewhat 
shorter than Eq. (155]) predicts, suggesting that the crust's thermal conductivity is higher than expected 



or that enhanced core neutrino cooling is taking place [94j. Another source, KS 1731-260 [95j, l96j, which 
produced a superburst lasting about 12 years, went into quiescence and has since shown a strong decline 
of X-ray luminosity, also suggesting a high thermal conductivity for the crust and rapid neutrino core 
cooling. Further work needs to be done, however, to test the assumption that the crusts of these sources 
are being heated to complete thermal equilibrium as otherwise shorter cooling times might be expected. 
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VI. COOLING NEUTRON STARS 

A. Observations of thermal emission 

The quantity inferred from thermal observations of a neutron star's surface is the so-called radiation 
radius 



#00 = R/\/l ~ 2GM/Rc 2 , (56) 

which results from a combination of flux and temperature measurements, both redshifted at the Earth 
from the neutron star's surface. For example, the un-redshifted Kirchoff's law for a black body, 

F =(f) 2 °T\ (57) 



becomes 



when redshifted to the Earth. Contours of the quantity R^ are displayed in Fig. [2j A measured value 
of Roo sets upper limits to both R and M, but without an independent estimate of mass or radius only 
limited constraints are possible. The major uncertainties involved in the determination of R^ include the 
distance, interstellar H absorption that is important near the peak and at lower energies of the spectrum, 
and details concerning the composition of the atmosphere and its magnetic field strength and structure. 
The most reliable measurements will originate from sources for which one or more of these uncertainties 



can be controlled. For example, the nearby sources RX J1856-3754 [97( and Geminga [98| have measured 
parallaxes, and the quiescent X-ray binaries in globular clusters for which reliable distances are (or will 
soon be) available. In addition, the X-ray binaries in globular clusters, having undergone recent episodes 
of accretion, are expected to have low magnetic field H-rich atmosperes, presumably, the simplest of all 
situations. 

A characteristic of these X-ray sources, in the cases in which the distances are small enough to allow 
detection of optical emission, is that the optical fluxes are a factor / ~ 5 — 7 times the amount expected 
from extrapolating the X-ray blackbody fits onto the Rayleigh- Jeans tail. This is one consequence of the 
neutron star atmosphere and its redistribution of the flux from a simple blackbody. The slope of the 
Rayleigh- Jeans tail is a measure of the optical temperature, and typically T opt ~ Tx/2. Therefore, in 
order to properly fit the overall flux distribution, a larger radius is needed than the X-ray fit alone would 
imply. The optical flux may be written as 

R Q pt \ 2 — . „ / Rj 



F ont oc47r ^ T, 



4nf(^f) 2 T x , (59) 



whereas the total stellar radius is a function of the effective optical and X-ray radii: 




R = jR 2 opt + R\ = R x Jl + -V ~ RxJl + 2f . (60) 



This explicitly shows that a factor / ~ 5 — 7 results in an increase in inferred radius of about a factor 
3-4 over that inferred from X-rays alone. Radii inferred from X-ray data alone are therefore suspect, and 
often much too small. 

Results from atmospheric fitting of the data for RX J1856-3754 and various globular cluster sources 
are displayed in Fig. [HJ An interesting feature of Fig. [8] is that if one accepts the results at face value, 
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FIG. 8: The same as Fig. [21 but displaying allowed mass and radius regions for thermally-emitting neutron stars. 
Top panel is for RX J1856-3754 (assuming a di stance of 120 pc) 971 ] . while the bottom panel shows results for 
globular cluster sources in M13 [99J], u> Cen [1QCH ] and 47 Tuc [101 ]. 
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only EOS curves that pass through all the permitted regions can be accepted. This would eliminate 
several relativistic field-theoretical EOSs such as GM3, FSU, MSO, MSI and MS2, while still permitting 
non-relativisitic potential models like AP3 and AP4. (It must be stressed, however, that the approach of 
Muller & Serot, of which the EOSs MSO, MSI and MS2 are typical examples, can be straightforwardly 
generalized to yield results that can accomadate the mass and radius constraints implied in Fig. 0) The 
Dirac-Brueckner relativisitic field-theoretical models like ENG and MPAl are allowed. Note that the 
surviving EOSs all support large masses, which could be crucial if the 2.1 M value for PSR J0751+1807 
remains robust. 



B. Implications for neutron star cooling 

For many years, it was common to adopt a standard cooling scenario, with relatively slow cooling, in 
which the dominant neutrino processes were the relatively slow modified Urea processes involving nucleons. 
However, neutrino production would become significantly more rapid if the symmetry energy and proton 
fraction is large or in the presence exo tic matte r such as hyperons, kaons or deconfined quarks, which 
would permit the direct Urea processes 



The direct Urea process involving baryons is the simple weak decays and the reverse reactions 

fli-> B 2 + £ + V e ; B 2 + l^ B 1 + u l , (61) 

where Bi and B 2 are baryons, and £ is a lepton, either an electron or a muon. When the decaying baryon 
is a neutron, the occurence of this process is contingent upon the presence of a sufficient number of protons 
(determined entirely by the density dependent symmetry energy) so that energy and momentum can be 
conserved simultaneously. 

In contrast to the direct Urea process, the modified Urea process includes a bystander baryon that 
facilitates the simultaneous conservation of energy and momentum. The additional particle slows the 
neutrino energy emission rate by a factor ~ (T/Tp) 2 , where Tp ~ 30 MeV is the Fermi energy at the 
equilibrium density n s . Since a typical temperature is of order 0.1 — 1.0 MeV, the modified Urea rate is 
millions of times slower than the direct Urea rate. However, the direct Urea process becomes possible, 
albeit with a smaller strength (matrix element), when a Bose condensate, deconfined quarks, or hyperons 
are present. 

The paradigm of slow standard (as opposed to rapid exotic) behavior stemmed from models of the 
nuclear force which had relatively small isospin dependences, such that the beta equilibrium proton and 
electron fractions would be very small and the direct Urea process among baryons would be forbidden. 
With the realization that newer versions of nuclear force models have greater isospin dependence and 
predict larger proton fra ction s has come the possibility that the direct Urea process can occur even in the 



absence of exotic matter [102| . Therefore, the detection of a rapidly cooling neutron star is not sufficient to 
unambiguously prove the existence of exotic matter in a neutron star interior. Compounding the difficulty 
of interpreting cooling observations is the role of superfluidity, which can turn off the direct Urea process 
near the superfluid's critical temperature. 

While it is outsi de the intended scope of this review to develop the details of neutron star cooling 



see, e.g., Ref. [l04 ] and references therein), it should be clear that it sensitively involves the internal 



composition. Inasmuch as the internal composition depends on the overall mass of the neutron star, 
an indirect measure of the mass of the star is possible. This is especially pertinent for t hose supernova 



remnants which have yet to yield detectable sources of thermal em issio n, such as 3C58 [105| , and supernova 
remnants G084.2-0.8, G093.3+6.9, G127. 1+0.5 and G315.4-2.3 [lOi}. If these remnants in fact contain 
neutron stars, the stars may have cooled rapidly due to one of the several direct Urea processes, which 
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is much more efficient than standard cooling due to the modified Urea process. Stars with high proton 
fractions may cool via the direct Urea process, and this depends upon the density dependence of the 
symmetry energy. A radiation radius determination for the same star would allow an extraction of the 
radius itself. 



VII. OTHER RADIUS OBSERVABLES 

Besides measurements of radiation radii from cooling neutron stars, ther limits to the radius could be 
set by pulse profiles, redshifts, explosions or accretion on neutron star surfaces with Eddington-limited 
fluxes, quasi-periodic oscillations (QPO's) from accreting neutron stars, and measurements of the moment 
of inertia of binary pulsars. We highlight these avenues in turn below. 



A. Pulse profiles 

Gravitational lensing strongly suppresses the amplitude of variations from pulsing (rotating) neutron 
stars. Essentially, gravitational light-bending allows an observer to see much more than the facing hemi- 
sphere. Therefore, analysis of pulsations in the thermal emission from pulsars allows limits to be placed 
on the ratio Mj R. The locally emitted flux is proportional the the local temperature to the fourth power, 
which varies since is controlled by the local magnetic field strength. Assuming dipole fields, it has been 
shown for many pulsars that pulse fractions comparable to the observed ones can be obtained only with 
stellar radii larger than those which are predicted by current models of neutron star struture, or with low 
stellar masses [59] . However, different field geometries can alter this conclusion. 

The same physics influences other pulsed emissions from neutron stars, such as Her X-l, which is an 
eclipsing X-ray emitting source with X-ray pulsations. The emission is observed to be occulted by the 
inner-disk edge which implies that the neutron star emits a pencil beam from its near pole and a fan 
beam from its far pole j6(J]. Gravitational light-bending strongly influences the emission, and Leahy [6l| 
derived M/R ~ 0.121 — 0.128 M Q km -1 , or z ~ 0.247 — 0.268, for this source. For the inferred mass of 



Her X-l, 1.29 - 1.59 M [34j, this translates into a radius range 10.1 < R/km < 13.1. This result may be 
model dependent, and more complex geometries and emissivities should be considered, but it nevertheless 
indicates the potential of the technique. 



B. Redshift 



Another observable is the redshift 



2GM\- 1/2 

1 H - 1 

Rc 2 J 



(62) 



Possibly the best case is that of the active X-ray burster EXO 0748-676 [107| for which a pair of resonance 
scattering lines consistent with Fe XXV and XXVI imply z ~ 0.345. Later, 45 Hz oscillations in the 
average power spect rum of 38 thermonuclear X-ray bursts were found and interpre ted as the neutron 

are consistent 
However, the 
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star spin frequency |l08j ]. Furthermore, the widths of the lines observed by Ref. 
with this spin frequency as lon g as the star's radius is between abou t 10 and 15 km 
spectral identification of Ref. 107] is not universally accepted |l09 |. Importantly, as long as the star's 
spin frequency is less than about 600 Hz, the measured z, which refers t o th e equatorial radius of the 
spinning star, is less than 2% different to that of the non-spinning star [lib] . More observations and 
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better modeling of the line profiles could lead to additional parameters involving M and R and serve as 
useful constraint. Other sources with lines have been claimed, but never confirmed with more sensitive 
instruments. 

Clearly, a simultaneous measurement of -Roo and z would determine both M and R: 

R = R^l + z)' 1 , 

M = ^ 00 (l + 2 )- 1 [l-(l + ^)- 2 ]. (63) 

The contour z = 0.35 is shown in comparison to typical mass-radius curves in Fig. [5J However, it is often 
the case that it is the quantity Roo/d that is determined, where d is the distance to the source, which is 
often uncertain. 

It is expected that planned X-ray missions, like Constellation-X which has a resolution and effective 
area several times larger than Chandra and XMM, will discover and precisely measure more spectral lines 
from neutron star surfaces. 



C. Eddington limit 

A major hurdle in determining radii from observations of thermal radiation from cooling neutron stars 
is that the distances, d, to the sources are generally not well known. Even if a source's redshift can be 
determined, an additional relation involving M, R and/or d is needed to independently find M and R. 
Burst sources, which are believed to involve thermonuclear X-ray bursts on neutron star surfaces, have 
long been known to have peak fluxes that are to order of magnitude comparable to the Eddington limited 
flux: 

cGM 



F edd ,oc = - 2GM/i?c 2 , (64) 

where k is the opacity. If the opacity is assumed to be dominated by electron scattering, one has 

k = 0.2(1 + X) cm 2 g" 1 , (65) 

where X is the hydrogen mass fraction. Following a burst, the source apparently relaxes to a quiescent 
state, which is believed to involve the radiation of thermal emission. Study of the X-ray spectrum can 
then yield a cooling flux and temperature Foo and observed at Earth: 

Foo/(oT^) = (Ttff/T^&dT^l - 2GM/Rc 2 )-\ (66) 
where T e ff is the effective temperature and is the observed color temp eratu re. The ratio T e ff/T 00 is a 



color correction factor that depends on the composition, gravity and T e ff [1111 . 

If the observed peak flux is approximately the Eddington flux, if the opacity and color correction 
factor are well-understood, if the quiescent flux observed in-between bursts is thermal emission with a 
measureable temperature T^, and if spectral features are observed which permit a determination of the 
object's redshift, these multiple observations then contain the information needed to identify uniquely the 
mass and the radius of a single star. To make this clear, combine Eqs. ([52"]) . (1641) and ([55]) : 

c 5 a 



R = [l-(l + z)- 2 ](l + z)- 3 
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c 3 \/a 



d = [l-(l + ^)- 2 ](l + ^)- 2 . (67) 

In the above, the factor a = ( T^/ Tpff^F^/jaTt) = (Roo/d) 2 . The burst source EXO 0748-676 is 



especially amenable to analysis [ll2 |. as this is one of the few sources with an estimated redshift (see 
above), z = 0.345. Its peak flux, to be identified with F ed d, is 2.25 ± 0.23 x 10~ 8 erg cm -2 s _1 , and the 
ratio FoajiaT 1 ^) is observed to be 1.14 ± 0.10 (km/ kpc) 2 . Using Eqs. ( 1671) and these observables, along 
with standard assumptions for k and T^/Teff, Ozel |l!2l | finds for this burst source M = 2.1 ± 0.28 M Q , 



R = 13.8 ± 1.8 km, and d = 9.2 ± 1.0 kpc. These errors do not include uncertainties in the assumed 
opacity or color-corrections. 



Note, however, that statistical studies [1131 ] indicate that the average maximum luminosity of X-ray 
bursts at infinity are about 60% larger than Eq. (1641) predicts. Various suggestions have been proposed 
to explain the super-Eddington luminosities in these bursts. Additionally, the value of k is as yet difficult 
to ascertain to good precision because the composition is uncertain and contributions from bound-free 
opacities are uncertain. Nevertheless, if these problems can be addressed, observations of this source could 
ultimately lead to the simultaneous determination of the mass and radius of a neutron star. In short, 
this is a promising avenue to pursue as it holds the promise to establish the mass and radius of the same 
neutron star. 

An interesting variation of the theme of the Eddington flux concerns the emissi on f rom neutron stars 



from accretion in low-mass X-ray binaries (LMXRBs). In the spreading layer model [1141 ]. accreting matter 
spreads from the equator to the poles on the neutron star surface, locally radiating at the Eddington flux. 
The total emission can approach the Eddington limit if the matter spreads all the way from the equator to 
the poles. Calculations indicate that the overal emitted spectrum is that of a diluted blackbody. Schemat- 
ically, the effective temperature of the radiation (i.e.,, the Eddington temperature Tem-, is determined by 
the balance between surface gravity and radiative acceleration: 



GM aT, 



4 

Edd 



RJl - 2GM/Rc 2 



k , (68) 



where n = 0.2(1 + X) cm 2 g 1 is the electron scattering opacity and X is the hydrogen mass fraction. The 
observed color temperature is 

T c = T Edd fJ 1/4 ^l-2GM/Rc 2 , (69) 

where i ~ 0.8 is the correction for centrifugal reduction of gravity and f c is a correction factor for spectral 
hardening. As a result, the neutron star radius can be found from 

2 £fy 2GM f 2GM\ 3 / 2 
R = 2^,— V - -R^J ■ (70) 

This is essentially equivalent to the middle formula in Eq. (IBTj) . In practice, modeling of the radiative 
transfer makes some corrections to the above, but it can be observed that the contours in an M — R 
diagram for fixed values of T c are nearly orthogonal to contours of Roo, and are less sensitive to variations 
in M. For 1.4 M neutron stars, and taking X ~ 1 and 2.3 < T c /keV < 2.5 (the observed range in 



LMXRBs), this method implies 12 < R/km < 15 [llSj. It is important to note that the spreading layer 



method only utilizes the color temperature of the source. If the same source could be observed after 
accretion ceases, the information obtained during accretion could be supplemented during quiescence to 
independently determine M and R as in the preceding analysis for EXO 0748-676. 
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D. Quasi-periodic oscillations from accreting stars 

Accreting matter in LMXRBs flows from the companion through its inner Lagrangian point into 
an accretion disc. Through viscous dissipation, the matter slowly loses angular momentum and spirals 
closer to the neutron star. Eventually, the matter reaches the inner edge of the accretion disc, which is 
thought to be close to the innermost circular stable orbit (ISCO), which in the Schwarzschild geometry 
is Risco — QGM/c 2 . Thereafter, the matter quickly spirals onto the neutron star's surface, emitting 
X-rays as g ravitational potential energy is converted into thermal energy. In some of these binaries, 25 
so far [ll6| . Fourier analysis of the X-ray emission shows quasi-periodic oscillations, or QPOs. Usually, a 
frequency of 200-400 Hz is observed, which is thought to be either the spin frequency of the star or half 
of that. In adddition, two high frequency kHz peaks are also seen, separated by the presumed neutron 
star spin frequency. Explanations o f the QPOs range from beat-frequency models connected wit h th e 
ISCO frequency and the star's spin [l!7 | to vertical epicyclic frequencies or resonant interactions [l!8 



Observations indicate that these QPOs rise slowly i n fre quency and suffer a drop in their quality factor 
as these frequencies saturate at their largest values [ligjj . The quality factor is the frequency divided by 



the peak's full width at half maximum. This behavior is what is expected if the drop is produced by the 
approach of accreting matter to the ISCO. It should be noted that all proposed models connect the upper 
kHz QPO with the ISCO orbital frequency to within a percent for neutron stars. 

Since the orbital frequency at the ISCO is a function of neutron star mass (if the star is rotating 
slowly), and the inner edge of the accretion disc must lie at or outside the ISCO, the inferred frequency 
thus provides a limiting mass for the neutron star. Inferred ISCO frequencies range from Vis cn ^ 1220 
Hz for 4U 1636-53, v ISCO ^ 1230 Hz for 4U 1608-52, and u ISCO ~ 1310 Hz for 4U 1728-34 [jig). The 



inferred upper limit to the masses, including the lowest order correction for neutron star spin, is [117 



/2200 Hz\ / cJ \ 

Ms b^)( 1+0 - 75 ^) Me - (71) 

where J = 2irfl is the angular momentum of a star spinning with frequency /. I is the moment of inertia 
(see §VII El) . Eq. (1711) is a slowly varying function of R, and the upper limit is shown by the nearly straight 
lines in the range 1.8-2.0 M Q in Fig. [6] for visco — 1300 Hz and neutron star spin frequencies / = 200 
and 400 Hz. The closer the observed QPO frequency is to that of the ISCO, the closer the neutron star 
mass must be to this upper limit. 



In addition, the neutron star radius must be smaller than the ISCO, which implies [117 



R 1950 Hz U 1+0 . 20 ^L )km . (72) 



visco ) V GM 

Radii limits for masses less than the upper limit Eq. (1721) scale with M 1//3 and are also shown in Fig. [61 
for the same visco an d /'s as for the mass limit. These radii limits are somewhat more restrictive than 
that implied by the 716 Hz pulsar J1748-2446ad. 

As upper limits to mass and radius, QPOs are not yet restrictive. However, if the measured frequencies 
are actually measurements of the ISCO, the implied masses are near 2 M Q . Interestingly, in the case of 
4U 1 636- 536, phase-resolved spectroscopy with the VLT suggest the neutron star is in the range of 1.6-1. 
M |l20j |. providing support for this interpretation and for the existence of high-mass neutron stars. 



E. Moments of inertia 



121 provides physicists with 



The discovery of the double-pulsar system PSR J0737-3039A & B |3J 
the most relativistic system yet. Perhaps most intriguingly, it could provide a measurement of spin-orbit 
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coupling that could eventually lead to a determination of the moment of inertia I a of star A 32| . Spin- 
orbit coupling advances the periastron of the orbit (apsidal motion) beyond the standard post-Newtonian 
advance as well as causi ng a precession of the orbital plane about the direction of the total angular 
momentum of the system 122| (geodetic precession). Since the masses of both stars are already accurately 
determined, and the moment of inertia can be expressed as a tight relation involving only M and R and 
no other EOS parameters, the measurement of I a could have enormous imp ortan ce. 



To d ate, t he best-determined moment of inertia is for the Crab pulsar [123 



This was based on an 

estimate |124j for the mass of the ionized portion of the Crab's remnant, 4.6 ± 1.8 M Q , which implies a 
lower limit to the Crab pulsar's moment of inertia of 97 ±38 M Q km 2 . This limit would rule out only very 
soft EOSs. 



According to Barker |l25| . the spin and orbital angular momenta evolve satisfying 
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G(Mi + 3M) 
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(73) 
(74) 



where the superscript SO refers to the spin-coupling contribution only (there are also first- and second- 
order post-Newtonian terms, 1PN and 2PN, respectively, unrelated to the spins, that contribute to this 
order). Here a is the semimajor axis of the effective one-body orbital problem (sum of the semi- major axes 
of the two stellar orbits), e is the eccentricity, M = Ma + Mb, and Mj refers to the mass of either binary 
component. To this order, one may employ the Newtonian relation for the orbital angular momentum: 
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where P is the orbital period and M = Ma + Mb- Then, from Eq. ( 1731) . the spin precession periods are 



P . 

Pi* 



2c 2 aPM{l - e 2 ) 
G(M — Mi) (Mi + 3M) ' 



(76) 



which are not identical for the two components unless they are of equal mass. Note that the spin precession 
periods are independent of the spins. Also note that if the spins are parallel to L, there is, first, no spin 
precession, and second, the spin-orbit contribution to the advance of the periastron is in the opposite sense 
to the direction of motion. 

are: 



The relevant observational parameters for PSR J0737-3039 extracted from Ref. [32 
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M B ~ 1.25 M , a/c = 2.93 s, e ~ 0.088 
P B ^ 2.77 s, P~ 0.102 day. 



(77) 



From these values, we observe that P p a — 74.9 yrs and P pB ~ 70.6 yrs. With these parameters, we can 
form the useful combinations 
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where Ia,so — ^a/(80 M km 2 ) is a typical value for the moment of inertia. 

The spin precession leads to changes in the directions of the pulsar beams. In many cases, this will 
lead to the periodic appearance and disappearance of the pulsar beam from the Earth. Second, since the 



33 



total angular momentum is conserved (to this order), the orbital plane will change orientation, which will 
be observed as a change in the inclination angle %. Damour & Schaefer 122j ] have considered the question 
of how these effects affect the timing of binary pulsars. For the change in i, 
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where Q{ is the angle between Si and L, and fa is the angle between the line of sight to pulsar i and the 
projection of Si on the orbital plane. The amplitude of the change in the orbital inclination angle i due 
to A's precession will be 
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This will cause a periodic departure from the expected time-of-arrival of pulses from pulsar A of amplitude 
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if one can assume that the orbital eccentricity is small. The facts that the orbit of PSR J0737-3039 is 
seen nearly edge-on and that S A is only slightly misaligned from L makes this a special case in which 
the amplitude of the timing change produced by the orbial plane precession will be extremely small, 
St a — (1.05 ± 1.33)7,4,80 A*s |S0|- Not only is the magnitude very small, but the large relative uncertainties 
in both cosi (| cos?| ~ 0.028±0.028) and in sin #,4 (sin #,4 ~ 0.22±0.17) combine to give a huge uncertainty 
in St a- We do not consider this effect further, but it could be relevant for future discoveries of binary 
pulsars with more favorable inclinations. 



For the advance of the periastron, the ratio of the spin-orbit and 1PN contributions is [122 
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In the case that \S A \ » \Sb\, on ly th e i = A term contributes substantially. For comparison, the ratio of 
the 2PN to 1PN contributions is [ill 
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where both % = A and i = B terms contribute. The second line of Eq. (183]) is valid in the case that 
Ma = Mb- Observational effects of the periastron advance are proportional to cosz and so the effect is 
maximized in this system. Furthermore, Eq. (182j) illustrates that uncertainties in the angle <pA are largely 
irrelevant in the case % ~ 90°. The 1PN periastron advance for PSR J0737-3039 is 
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or 0.294 radians per year. The periastron advance ratio works out to be [80 
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FIG. 9: The moment of inertia scaled by M 3 / 2 as a function of stellar mass M for EOSs described in The 
shaded band illustrates a ±10% error on a hypothetical //M 3//2 measurement with centroid 50 km 2 M Q 1//2 ; the 
error bar shows the specific case in w hich the mass is 1.34 M Q with essentially no error. The dashed curve labelled 



"Crab" is the lower limit derived by [1231 ] for the Crab pulsar 



However, a practical difficulty remains. The spin-orbit contribution depends upon the individual masses, 
not just the total mass which is known to extremely high precision, requ iring; for our purposes the deter- 
mination of essentially three post-Newtonian parameters, as Damour |l22| has emphasized, to an accuracy 



of a part in 10 5 . The accuracy of aesini at present is a few parts in 10 3 , and the determination of I a 
depends on its refinement by a factor of 100 or more. There is reason to be optimistic that this can occur 
within a decade (80| . 

A sufficiently accurate measurement of / say, to within 10%, would usefully discriminate among families 
of EOSs because the mass will be known to great precision. This is displayed in Fig. [9j in which the 
quantity J/M 3 / 2 (in order to reduce the vertical range) is plotted versus M. The shaded band illustrates 
a hypothetical measured value with error bars of ±10%. Only curves that pass through the line with the 
error bar would be allowed. In this hypothetical case, the vast majority of the E OSs illustrated would not 



be adequate. For comparison, the lower limit derived for the Crab pulsar [1231 ] . which rules out only the 
softest EOSs, is also shown in Fig. [91 This limit is not nearly so restrictive. 

The fit in Eq. ([35]) for EOSs previously cited as supporting masses greater than about 1.6 M Q is 
reasonably tight. This relation may be inverted to yield a radius estimate given I and M. Furthermore, 
since / oc R 2 , the inferred relative error of R using this procedure would be approximately half the relative 
measured error of I itself. 
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VIII. NEUTRINO EMISSION FROM A PROTO-NEUTRON STAR 



A proto-neutron star (PNS) is born in the aftermath of a successful supernova explosion resulting 
from the gravitational collapse of the core of a massive star. During the first tens of seconds of evolu tion, 
nearly all (~ 99% ) of the remnant's binding energy is radiated away in neutrinos of all flavors [si, Il2(j 
I27L EH, E3, 130 ]. The neutrino luminosities and the emission timescale are controlled by several factors 
including the total mass of the PNS and the opacity of neutrinos at supranuclear density, which depend 
on the star's composition and the EOS of strongly interacting dense matter. One of the chief objectives in 
modeling PNSs is to infer their internal compositions through future su perno va neutrino signals detected 



in neutrino observatories like SuperK, SNO and others, including UNO [131 



A. The birth of proto-neutron stars 



The evolution of a PNS proceeds through several distinct stages and with various outcomes 1 133 



13J], as shown schematically in Fig. [TU1 Immediately following core bounce and the passage of a shock 
through the PNSs outer mantle, the star contains an unshocked, low entropy core of mass ~ 0.7 M Q in 
which neutrinos are trapped (stage I in the figure). The core is surrounded by a low density, high entropy 
(5 < s < 10) mantle that is both accreting matter from the outer iron core falling through the shock and 
also rapidly losing energy due to electron captures and thermal neutrino emission. The mantle extends up 
to the shock, which is temporarily stalled about 200 km from the center prior to an eventual explosion. 

After a few seconds (stage II), accretion becomes less important if the supernova is successful and the 
shock has ejected the stellar envelope. Extensive neutrino losses and deleptonization will have led to a loss 
of lepton pressure and the collapse of the mantle. If enough accretion has occurred, however, the star's 
mass could increase beyond the maximum mass capable of being supported by the hot, lepton-rich matter. 
If this occ urs, the remnant collapses to form a black hole and its neutrino emission is believed to quickly 



cease 
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Neutrino diffusion deleptonizes the core on time scales of 10-15 s (stage III). Diffusion time scales 
are proportional to -R 2 (cA„) -1 , where R is the star's radius and X u is the effective neutrino mean free 
path. This generic relation illustrates how both the EOS and the composition influence evolutionary time 
scales. The diffusion of high-energy (200-300 MeV) z/s from the core to the surface where they escape 
as low-energy (10-20 MeV) z/s generates heat (a process akin to joule heating) through neutrino-matter 
interactions. The core's entropy approximately doubles, producing temperatures in the range of 30-60 
MeV during this time, even as neutrinos continue to be prodiguously emitted from the star's effective 
surface, or neutrinosphere. 

Strange matter, in the form of hyperons, a Bose condensate, or quark matter, suppressed when neutri- 
nos are trapped, could appear at the end of the deleptonization. Its appearance would lead to a decrease in 
the maximum mass that the internal pressure of matter is capable of supporting a gain s t gra vity, implying 
metastability of the neutron star and another chance for black hole formation [1331 . Il34j . This would 



occur if the PNSs mass, which must be less than the maximum mass of hot, lepton-rich matter (or else 
a black hole would already have formed), is greater than the maximum mass of hot, lepton-poor matter. 
However, if strangeness does not appear, the maximum mass instead increases during deleptonization and 
the appearance of a black hole would be unlikely unless accretion in this stage remains significant. 

The PNS is now lepton-poor, but it is still hot. While the star has vanishing net neutrino number, 
thermally produced neutrino pairs of all flavors dominate the emission. The average neutrino energy slowly 
decreases, and the neutrino mean free path increases. After approximately 50 seconds (stage IV), A ~ R, 
and the star finally becomes transparent to neutrinos. The neutrino luminosity rapidly decreases beyond 
this time. Since the threshold density for the appearance of strange matter decreases with decreasing 
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FIG. 10: The main stages of evolution of a neutron star, from Ref. 133 ]. Shading indicates approximate relative 
temperatures. 

temperature, a delayed collapse to a black hole is still possible during this epoch. 

Following the onset of neutrino transparency, the core continues to cool by neutrino emission, but the 
star's crust remains warm and cools less quickly. The crust serves as an insulating blanket which prevents 
the star from coming to complete thermal equilibrium and keeps the surface relatively warm (T m 3 x 10 6 
K) for up to 100 years (stage V). The temperature of the surface after the interior of the star becomes 
isothermal (stage VI) is determined by the rate of neutrino emission in the star's core and the composition 
of the surface. 



B. Theoretical expectations 



To understand what aspects of the EOS and str ucture can be probed by neutrinos, we examine 
some analytic models for proto-neutron star evolution |l33l 136 ] . For clarity and simplicity, we employ 
Newtonian gravitation, as this does not affect the qualitative conclusions we will draw. We will assume 
that the neutrino distribution function is well-approximated by a Fermi-Dirac distribution, so the neutrino 
number density is n u = J °° n u (E u )dE u , where 
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(86) 
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We also make the diffusion approximation, in which both number and energy fluxes are driven by density 
gradients: 
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(87) 



where F v and L v are the net electron neutrino number flux and total neutrino luminosity, respectively, \ v 
(Ap)and A^ are the mean free paths for number and energy transport, respectively, and 6j is the neutrino 
energy density i = [u e , u e , P M , 

number Yl = Y e + Y v — Y p and the total entropy per baryon are 

0Y L 1 



z/ e , z/ e , v„, z/„, i/ r , z/,.] refers to the neutrino species. The net changes in electron lepton 
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where j = n,p, e, v e . 

The net neutrino mean free paths \ v and \ l E are due to both absorption and scattering and should 
also include inverse processes. In what follows, we will want to have an approximation for X u under 
degenerate conditions and when electron neutrino absorption dominates. Similarly, we will want to have 
an approximation for \ % E under non-degenerate conditions and when it is dominated by scattering. In 
that case, the transport is dominated by all neutrino types other than u e , and one can replace the sum 
in Eq. (188|) by a single term involving the scattering mean free path. Ref. 136l | shows that suitable 
approximations are then 
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where Ao — 50 cm, Ac ~ .2 km and T ^ 
function under degenerate conditions: 
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10 MeV. Finally we note a property of the Fermi distribution 
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It is convenient to examine two distinct periods of proto-neutron star evolution: the deleptonization 
era and the cooling era. During the first of these, the transport within the star is dominated by degenerate 
electron neutrinos propagating through degenerate matter. Only near the surface do these conditions 
break down. We can use the first of Eqs. ( IHHi) to examine the loss of leptons from the star. Using Eq. 
(EHD for \ v , Eq. (EZD for F v , and Eq. USD), one has 



OnY L 
Ot 



cA 



1 



]Tr 2 (hc) 3 r 2 Or 



T \ 2 0f4 
Or 



-) 



(91) 



Using 0Y L /0Y u ~ (0Y L /0Y u ) ~ 5 and 6ii 2 nY u (hc) 3 = for degenerate neutrinos, ignoring spatial 
variations of T since the core is nearly isothermal, and separating the spatial and temporal variations 
using fil = {il ip(x)(j)(t), this can be rewritten as 
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where 



x = r H> rD = ^^J [tJ {wj a - (93) 

T c is the assumed uniform temperature. The quantity X\ is the value of x for which i/j first vanishes. 
The second equality of Eq. ( l92l) is the Lane-Emden equation of index 1 for which the solution is ip(x) = 
sin(x)/x and has the outer boundary x\ = ir. The first equality of Eq. ( l9Ti) is an exponential decay 
<j)(t) = exp(— £/td) and td is clearly the diffusion timescale. It is interesting that td is independent of 

fi U fi but does depend sensitively on the ambient temperature T c . Since one expects T c oc n^ nteT at fixed 
entropy, td will scale with the mass of the remnant like M 4 / 3 since R changes very slowly with M. 
From Eq. flH7|) . one can determine the neutrino flux 

At the outer boundary, where xdip/dx = — 1, the emerging flux is then 

where T e (t) is thereby defined as the effective blackbody temperature, and F t is the standard Fermi integral 
of index i. We assume that the effective value of the chemical potential when the neutrinos decouple from 
the matter is fi u<e < T. The effective temperature is therefore 

Te(t)^(j^) 1/3 (^) 2/3 ^(t) 1/3 - (96) 
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The average emergent neutrino energy is then 

< E v >~ T e (t)(F 3 (0)/F 2 (0)) ~ 3.15T e (t) , (97) 

and the star's emergent luminosity will be At[R 2 F v < E v >. The effective value of the chemical potential 
at decoupling can be estimated using the Eddington approximation, which gives 

* - -I'i (98) 

where A e is the neutrino mean free path at the decoupling point. From this, one can show that fi e ^ T e so 
we are justified in using zero arguments in the above Fermi integrals. 

During the beginning of the proto-neutron star phase, R ~ 20 km and T c ~ 20 MeV. One thus obtains 
r D = 16 s, F v (R,0) = 2.5 x 10 42 neutrinos cm" 2 s -1 , T e ~ 2.9 MeV and < E v >~ 9.2 MeV. The initial 
electron neutrino luminosity due to deleptonization is L u (R,0) = 2.1 bethe s -1 . We remark that the 
initial luminosity burst predicted by detailed models is much larger, around 100 bethe s _1 , but these are 
neutrinos originating from the outer mantle of the proto-neutron star near the neutrinosphere. Within 
a few seconds, as the mantle deleptonizes and collapses, the emergent luminosity becomes dominated by 
core emission. 

Energy transport is also important during the deleptonization stage. Noting that 

HjdYj = (-/i n +Hp + n e - [i v )dY e + [i v dY L , (99) 
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the first term on the right-hand side vanishes in beta equilibrium. For electron neutrinos, which dominate 
in the interior, the mean free paths are so small that the first term on the right-hand side of Eq. ([55]) is 
negligible during deleptonization compared to the second term. We therefore have 
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This shows that continuous replenishment of diffusing neutrinos enhances the growth of entropy in the 
star. For degenerate matter, s ~ aT where a ~ 0.05 MeV -1 is the supra-nuclear specific heat. Integration 
leads to 
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where Sf and s, are the final and initial entropies. For Si ~ 1 one obtains s/ ~ 2.5. The final entropy 
scales roughly as M 1 / 6 . Thus it is clear that the central portions of the star heat significantly during 
deleptonization. 

Following deleptonization, when the central temperature reaches a maximum, core cooling begins. By 
this time, fi u < T and its gradient is also small. The second term in the second of Eq. ( |88|) becomes 
negligible. Assuming that the energy transport mean free paths of all flavors of neutrinos are roughly equal, 
we can replace the sum in the second of Eq. (I57|) by a factor 6, and the neutrino luminosity becomes 
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where we used the fact that J °° EfdE = ir 2 T 2 /12. During the cooling phase, the density gradients in 
the star's center are small, and one can remove the density dependence from the integral in the energy 
transport equation, which beomes 
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We let be the density in the star's center. The matter specific heat is dominated by baryons, so that 

(104) 
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where a ~ 0.1 is the level density of nucleons at n s . Once again, we can seek a separable solution with 

T = T^(x)(f)(t): 
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where x\ ~ 19 is the eigenvalue for the index 2 Lane-Emden equation. The temporal solution is trivial, 
0=1 — t/rc- T* is the core temperature at the end of deleptonization and the beginning of cooling (t = 
in the immediate above). 

We can establish the emergent luminosity and define an effective temperature with 
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FIG. 11: The binding energy per unit mass of various EOSs as a function of M/R. The lighter curves are EOSs 
for nucleonic stars as well as those containing hyperons, Bose condensates and deconfined quark matter (see Ref. 
0] for details). The five heavier curves are for analytic solutions of GR structure equations. The shaded band is 
the fit with error extents as given in Eq. (|108p . Figure taken from Ref. [fjj]. 



For the index 2 Lane-Emden equation, —(xdip/dx) Xl ~ 0.56. Inserting values for T* = 50 MeV and 
n* = 4n s , we find r c ~ 18 s, L U (R, 0) ~ 11 bethe s _1 , T e (0) ~ 5.5 MeV, and < E v >~ 17 MeV. Note that 
T e at the beginning of the cooling era is nearly twice as large as at the beginning of deleptonization. This 
is attributable to core heating which occurs during deleptonization and results in a steady increase in the 
average energy of emitted neutrinos during the deleptonization era. Finally, note that like 
Also T e oc Tl' 4 oc M 1 / 6 . 

In summary, there are a multitude of observables that can provide details of the proto-neutron star: 
timescales, average energies, and neutrino light curves. Of course, the unknown details of the proto-neutron 
star such as its mass, EOS and opacities are not related to these variables in unambiguous ways. However, 
one observable that does not depend upon opacities can be considered: the total neutrino emission which 
equals the binding energy released in the gravitational collapse to form the neutron star. In Newtonian 
gravity, the binding energy of a uniform sphere is BE = 3GM 2 /(5R). Ref. determined an approximate 
fit to the binding energies for EOSs that is valid for EOSs which permit maximum masses larger than 
about 1.65 M : 

BE/M ~ (0.60 ± 0.05)/3(l - /5/2)" 1 , (108) 
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where (3 = GM/Rc 2 . The binding energies of neutron stars formed from a variety of EOSs is shown, 
together with the above fit, in Fig. [TTJ Additionally are shown binding energies for five analytic solutions 
of GR applicable to neutron or self-bound stars 0. With the detection of thousands or millions of 
neutrinos, the binding energy could be determined to better than a percent, and the tightness of this fit 
would imply a correspondingly tight relation between M and R to be established. 

C. Model calculations 
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FIG. 12: The evolution of the average energy and total luminosity of neutrino s in P NSs composed of baryons only 
(left panel) and baryons and hyperons (right panel). The figure is from Ref. 
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Neutrino signals from PNSs depend on many stellar properties, including the mass; initial entropy, 
lepton fraction and density profiles; and neutrino opacities. In Figs. [12] -[HI the de pendence of ne utrino 
emission on PNS characteristics are shown from the detailed study of Pons et al. |l28l . 129 . 130| . The 



generic results (see Fig. [T2|) are that both L v and < E u > increase with increasing mass (il 128 . < E u > 



for all flavors increases during the first 2-5 seconds of evolution, and then decreases nearly linearly with 
time. For times larger than about 10 seconds, and prior to the occurrence of neutrino transparency, the 
L v decays exponentially with a time constant that is sensitive to the high-density properties of matter. 
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FIG. 13: The evolution of the total neutrino luminosity for npQ PNSs. Shaded bands illustrate the limiting 
luminosities corresponding to a count rate of 0.2 Hz, assuming a supernova distance of 50 kpc for 1MB and 
Kamioka, and 8.5 kpc for SNO and SuperK. The widths of the shaded regions represent uncertainties i n th e 
average neutrino energy from the use of a diffusion scheme for neutrino transport. Figure taken from Ref. 13ol ] . 



Significant variations in neutrino emission occur beyond 10 seconds: L v is larger during this time for stars 
with smaller radii and with the inclusion of hyperons in the matter. Finally, signi fican t regions of the stars 
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become convectively unstable during the evolution, as several works have found 

The main effect of the larger mean free paths produced by in-medium corrections 127 . 138| is that the 
inner core deleptonizes more quickly. In turn, the maxima in central temperature and entropy are reached 
on shorter timescales. In addition, the faster increase in thermal pressure in the core slows the compression 
associated with the deleptonization stage, although after 10 s the net compressions of all models converge. 
The relatively large, early, changes in the central thermodynamic variables do not, however, translate into 
similarly large effects on observables such as L v and < E v >, relative to the baseline simulation. It is 
especially important that at and below nuclear density, the corrections due to correlations are relatively 
small. Since information from the inner core is transmitted only by the neutrinos, the time scale to 
propagate any high density effect to the neutrinosphere is the neutrino diffusion time scale. Since the 
neutrinosphere is at a density approximately 0.01n s , and large correlation corrections occur only above 
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n s /3 where nuclei disappear, correlation corrections have an effect at the neutrinosphere only after 1.5 
s. However, the corrections are still very important during the longer-term cool i ng st age, and result in a 
more rapid onset of neutrino transparency compared to the Hartree results 127t 138l . 




FIG. 14: Lifetimes of metastable stars versus the PNS baryon mass Mb- Thick lines denote cases in which 
the maximum gravitational masses of cold, catalyzed stars are near 1.45 M , which minimizes the metastability 
lifetimes. The thin lines for the npQ and npH cases are fo r EOSs with larger maximum gravitational masses (1.85 



and 1.55 M , respectively.) Figure taken from Ref. 13(J ] 



A comparison of the signals observable with different detectors is shown in Fig. [131 which displays L u 
as a function of baryon mass Mb for stars containing quarks in their cores. In the absence of accretion, Mb 
remains constant during the evolution, while the gravitational mass Mq decreases. The two upper shaded 
bands correspond to estimated SN 1987A (50 kpc distance) detection limits with KII and 1MB, and the 
lower bands correspond to estimated detection limits in SNO, SuperK, and UNO, for a Gala ctic supernova 
(8.5 kpc distance). The detection limits have been set to a count rate dN/dt = 0.2 Hz 129[. It is possible 
that this limit is too conservative and could be lowered with identifiable backgrounds and knowledge of 
the direction of the sign al. The width of the bands represents the uncertainty in < Ep e > due to the 



diffusion approximation [128l . Il29l . Il30| . It appears possible to distinguish between stable and metastable 



stars, since the luminosities when metastability is reached are always above conservative detection limits. 



D. Metastable proto-neutron stars 



Proto-neutron stars in whi c h str angeness appears following deleptonization can be metastable if their 



masses are large enough fl39l . [lib] . This stems from the interesting behavior that the expression of 



strangeness is suppressed in hot, lepton-rich stars because of the large value of /i e . However, as stars 
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deleptonize, /i e falls and it is possible that strangeness eventually appears in the star's core. Generally, 
this softens the equation of state, causing a decrease in the effective maximum mass. This leads to a 
possibly important diagnostic that could shed light on the internal composition of neutron stars: the 
abrupt cessation of the neutrino signal when a deleptonizing strangeness-containing star suddenly finds 
that its mass exceeds the maximum mass. This would be in contrast to a normal star of similar mass for 
which the signal continues to fall until it is obscured by the background. In Fig. [141 the lifetimes for stars 
containing hyperons (npH), kaons (npK) and quarks (npQ) are compared |l29l |. In all cases, the larger 
the mass, the shorter the lifetime. For the kaon and quark PNSs, however, the collapse is delayed until 
the final stage of the Kelvin- Helmholtz epoch, while this is not necessarily the case for hyperon-rich stars. 
In addition, there is a much stronger mass dependence of the lifetimes for the hyperon case. 

Clearly, the observation of a single case of metastability, and the determination of the metastability 
time alone, will not necessarily permit one to distinguish among the various possibilities. Only if the 
metastability time is less than 10-15 s, could one decide on this basis that the star's composition wa s tha t 



of npH matter. However, as in the case of SN 1987A, independent estimates of Mb might be available [141 



In addition, the observation of two or more metastable neutron stars might permit one to differentiate 
among these models. This highlights the need for breakthroughs in lattice simulations of QCD at finite 
baryon density in order to unambiguously determine the EOS of high density matter. In the meantime, 
intriguing possible extensions of supernova and PNS simulations with np Q an d npK matter include the 
consideration of heterogenoeus structures and quark matter superfluidity |l42l . 



IX. MAGNETIC MATTERS 



Recent investigations of the effects of ultra-strong magnetic fields (B > 10 14 Gauss) on neutron stars 
are spurred by several independent arguments that link the class of soft 7— ray repeaters and perhaps 
certain ano malous X-ray puls ars with neutron stars having ultra strong magnetic fields - the so-called 



magnet ars (l43l \ul \uH H 



In addition, some soft 7— ray repeaters directly imply, from their periods 
and spin-down rates, surface fields in the range 2 — 8 x 10 14 G auss. (See Table I from Cardall et. al. 



1471 ] for a summary of observations.) Kouveliotou et al. |148j argue from the population statistics of 



soft 7— ray repeaters that magnetars constitute about 10% of the neutron star population. While some 
observed white dwarfs have large enough fields to give ultra-strong neutron star magnetic fields through 
flux conservaton, it does not appear likely that such isolated examples could account for a significant 
fraction of ultra-strong field neutron stars. Therefore, an alternative mechanis m seems n ecessary for the 
creation ultra-strong magnetic fields in neutron stars. Duncan & Thompson [l49 . 150 ] suggested that 



large fields (up to 3 x 10 x (1 ms/P;) Gauss, where Pi is the initial rotation period) can be generated in 
nascent neutron stars through the smoothing of differential rotation and convection. 

The effects of magnetic fields on the E OS a t low densities , relevant for neutron star crusts, has been 
extensively studied (see for example, Refs. 151 . 152 . 153 . 154j ). Only a handful of w orks have consi dered 



the effects of very large magnetic fields on the EOS of dense neutron star matter [155J, Il56l . Il57l |. Lai 



and Shapiro [15 XI ] considered non-interacting, charge ne utral, beta-equilibrated matter at subsaturation 
densities, whereas Chakrabarty and co-authors 1551 . Il56| studied dense matter including interactions using 
a field-theoretical description. These authors found large compositional changes in matter induced by 
ultra-strong magnetic fields due to the quantization of orbital motion. Acting in concert with the nuclear 
symmetry energy, Landau quantization substantially increases the concentration of protons compared to 
the field-free case, which in turn leads to a softening of the EOS. This lowers the maximum mass relative 
to the field-free value. In these works, however, the electromagnetic field energy density and pressure, 
which tend to stiffen the EOS, were not included. In addition, changes i n the general relativistic structure 



induced by the magnetic fields (studied in detail by Bocquet et al. [I58j who, however, omitted the 
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compositional changes in the EOS due to Landau quantization) were also ignored. Calculations including 
the combined effects of the magnetic fields on t he EOS and on the general relativistic structure have been 
reported by Cardall, Prakash & Lattimer 147 . 



The most intriguing questions are: 

1. What is the largest frozen-in magnetic field a stationary neutron star can sustain?, and, 

2. What is the effect of such ultra-strong magnetic fields on the maximum neutron star mass? 



A. Magnetic effects on the EOS 

The answers to both of these questions hinge upon the effects strong magnetic fields have both on 
the EOS of neutron-star matter and on the structure of neutron stars. The magnitude of the magnetic 
field strength B nee ded to dramatically affect neutron star structure directly can be estimated with a 



dimensional analysis [1511 ] equating the magnetic field energy E b ~ (47iR 3 /3)(B 2 /8ir) with the gravitational 



binding energy Eb.e. ~ 3GM 2 /5R, yielding the so-called virial limit 

B ~ 1.4 x io 18 ( -4^r) (— t— ) 2 Gauss ' ( 109 ) 



where M and R are, respectively, the neutron star mass and radius. 

The magnitude of B required to directly influence the EOS can be estimated by considering its effects 
on charged particles. Charge neutral, beta-equilibrated, nucleonic matter contains both negatively charged 
leptons (electrons and muons) and positively charged protons. Magnetic fields quantize the orbital motion 
(Landau quantization) of these charged particles. Relativistic effects become important when the particle's 
cyclotron energy eh B/(mc) is comparable to it's rest mass energy. The magnitudes of the so-called critical 
fields (although there is nothing critical about this) are 

B e c = (hc/e) A e ~ 2 = 4.414 x 10 13 Gauss , 
B» = (m^/m e ) 2 B e c = 1.755 x 10 18 Gauss, 

B v c = (m p /m e ) 2 B e c = 1.487 x 10 20 Gauss, (110) 

for the electron, muon and proton, respectively (A e = h/m e c ~ 386 fm is the Compton wavelength of 
the electron). When the Fermi energy of the proton becomes significantly affected by the magnetic field, 
the composition of matter in beta equilibrum is significantly affected. In turn, the pressure of matter is 
significantly affected. In neutron star matter this situation occurs when B* = BjB c e ~ 10 5 , and will lead 
to a general softening of the EOS as shown by Broderick, Prakash & Lattimer 159| . 



In neutron stars, magnetic fields may well vary in strength from the core to the surface. The scale 
lengths of such variations are, however, usually much larger than the microscopic magnetic scale l m , which 
depends on the magnitude of B. For low fields, for which the quasi-classical approximation holds, 

l m ~ (A e 2 /B*)(37r 2 n e ) 1/3 « 10 5 (n e K) 1/3 /£* fm > ( m ) 

where n e is the number density of electrons and n s is the normal nuclear saturation density (about 0.16 
fm -3 ). For high fields, when only a few Landau levels are occupied, 

l m ~ 2ix 2 n e (X 2 jB*) 2 « 7 x 10 9 (n e /n s )/B* 2 fm. (112) 

In either case, the requirement that R » l m is amply satisfied; hence, the magnetic field B may be 
assumed to be locally constant and uniform as far as effects on the EOS are concerned. 
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In non-magnetic neutron stars, the pressure of matter ranges from 2 — 5 MeV fm 3 at nuclear density 
to 200 - 600 MeV fm~ 3 at the central density of the maximum mass configuration, depending on the EOS 



133l | . These values may be contrasted with the energy density and pressure from the electromagnetic field: 



B 2 

P f = — = 4.814 x 10' 8 B* 2 MeV fm 



113) 



Note that e = P f = 1 MeV fnT' 3 for a B* = 4.56 x 10 s or B ~ 2 x 10 14 Gauss. Thus, the field 
contributions can dominate the matter pressure for B* > 10 4 at nuclear densities and for B* > 10 5 at 
the central densities of neutron stars, and must therefore be included whenever the field dramatically 
influences the star's composition and matter pressure. 

In strong magnetic fields, contributions from the anomalous magnetic moments of the nucleons must 
also be considered. Experimentally, 



1| for the proton and K n — A*iV^ for the neutron . 



;ii4) 



where fi^ is the nuclear magneton and g p = 5.58 and g, 
and neutron, respectively. The energy 



: —3.82 are the Lande g-factors for the proton 
Kn + Kp\B ~ 1.67 x 10~ 5 B* MeV (115) 



measures the changes in the beta equilibrium condition and to the baryon Fermi energies. Since the Fermi 
energies range from a few MeV to tens of MeV for the densities of interest, the contributions from the 
anomalous magnetic moments become significant for B* > 10 5 . In fact, complete spin polarization of the 
neutrons occurs when 



\B 



(67r 2 n r 



X2/3 



4m T 



lie) 



At nuclear density, this leads to B* ~ 1.6 x 10 5 or B ~ 7.1 x 10 18 Gauss. Such spin polarization results in 
an overall stiffening of the EOS (due to the in creas ed degeneracy pressure of neutrons) that overwhelms 
the softening induced by Landau quantization 159l . 

Inasmuch as the anomalous magnetic moments of stra ngen ess-bearing hyperons are similar in magni- 
tude to those of nucleons (for a compilation, see Table 1 of 160], magnetic field interactions will be present 
if hyperons appear in dense matter. The influence of magnetic fields on the EOS of matter containing 
strangeness-bearing hyperons is characterized by a suppression of hyperons relative to the field-free case 
as both Landau quantization and the anomalous magn etic moment interactions serve to postpone the 



densities at which hyperons appear in dense matter [160 



AuM. and A4 is the magnetization [161 
B 2 , are introduced into the structure equations 147 



In magnetized matter, the stre ss en ergy tensor contains terms proportional to H B, where H = B + 

Thus, extra terms, in addition to the usual ones proportional to 
The mag netiz ation in a single component electron 
for neutron star crust matter. The 
generalization of this formulation to the case of interacting multicomponent matter with and without the 
effects of the anomalous magn etic moments leads to the result that deviations of H from B occur only for 



gas has been studied extensively by Blandford and Hernquist [162 



field strengths B* > 10 5 [159 



In a magnetic field, the energy of an electron receives contributions from vacuu m pol arization effects 
as well. Up to one loop level, these contributions have been calculated by Schwinger 163| with the result 
E% ac = m e c 2 n e , where the dimensionless factor K e (arising from the anomalous magnetic moment of the 
electron) takes on sign changing values in the weak and strong field limits: 



H 



4tt 



weak field result 
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K* e f = — (hi2B*) , strong field result , (117) 

where a e ~ 1/137 is the electromagnetic fine structure constant. One can set \k™*\ = 1 and k s J = 1 in 
order to gauge the fields required to induce a change in the electron Fermi energy (or chemical potential 
at zero temperature) by an amount m e c 2 . The results are: 

Air 

B* = 1.7 x 10 3 or B ~ 10 17 Gauss , for weak fields 

Q/ 1 ■ 

= — e^ 2n/ae ~ 2.8 x 10 12 or B ~ 10 26 Gauss , for strong fields . (118) 
2n 

Insofar as the effects of symmetry energy (with shifts of tens of MeV's), Landau quantization and the 
anomalous magnetic moments of baryons overwhelm the influence on proton fractions, vacuum polariztion 
effects of the electron are not significant except for extremely large fields. 

In addition, it must be borne in mind that for super-strong fields exceeding the proton critical field 
of ~ 1.5 x 10 20 Gauss, the energy density in the field would be significantly larger than the proton's mass 
energy density. At the proton critical field, the field energy density B 2 /(8tv) ~ 0.5 GeV fm -3 is comparable 
to the energy density of the proton ~ 1 GeV fm -3 , which calls for taking the compositeness of the proton 
seriously. 

Figure [15] shows typical results when the effects of Landau quantization and anomalous magnetic 
moments are included in nucleonic matter. The main lesson to be learned from these results is the 
softening of the EOS due to Landau quantization, which is, however, overwhelmed by stiffening due to 
the incorporation of the anomalous magnetic moments of the nucleons. These effects become significant 
only for fields in excess of B* ~ 10 5 , for which neutrons become completely spin polarized. Note that this 
field strength is substantially less than the proton critical field. In addition, the inclusion of ultra-strong 
magnetic fields leads to a reduction in the electron c hemi cal potential and an increase in proton fraction. 
Similar results obtain upon the inclusion of hyperons |l60j |. These compositional changes have i mplic ations 
for neutrino emission via the direct Urea process and, thus, for the cooling of neutron stars 156|. The 
magnetization of the matter never appears to become very large, as the value of \H/B\ never deviates 
from unity by more than a few percent. 



B. Magnetic effects on neutron star structure 

We turn now to address whether or not stable stellar configurations can exist in which the magnetic 
field is large enough (B > 5 x 10 18 G) that the properties of m atter are significantly affected by the magnetic 



field. One step in this direction was undertaken in Ref . [147] , in which the limits of hydrostatic equilibrium 
for axially-symmetric magnetic fields in general relativistic configuratio ns w e re a nalyzed, including the 
effects of the magnetic field on the EOS. As discussed in detail in Ref. 147 . 158j |. in axially symmetric 



field configurations with a constant current function, the magnetic field contributes a centrifugal-like 
contribution to the total stress tensor. This can be understood by noting that for the field geometry 
considered (namely axisymmetric) , a superconducting fluid can move along field lines but not across 
them. Thus the "pressure" associated with the magnetic field will only act equatorially and not vertically. 
This flattens an otherwise spherical star, and for large enough fields, decreases the central (energy) density 
even as the mass is increased. For large enough fields, the star's stability is eventually compromised. 

The question of when is the magnetic field too large is answered to an extent by the results shown in 
Fig. [H3 Results shown here are for the GM3 EOS, both excluding (lower panel) and including (upper 
panel) hyperons. As with rotation, magnetic fields allow neutron stars with a particular EOS and baryon 
number to have larger masses and equatorial radii compared to the field-free case. The maximum mass 
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FIG. 15: Matter pressure -P m , nucleoli Dirac effective mass vri^j vn n ^ and concentrations Y% — nijnij as functions of 
the density u = rib/n s (left panels; n s = 0.16 fm -3 is the fiducial nuclear saturation density) and magnetic field 
strength B* = BjB c e (right panels; B c e = 4.414 x 10 13 Gauss is the electron critical field), for the model GM3. The 
inset in the upper left panel shows P m as a function of the matter energy density e m . The curve labeled P* in the 
upper right panel shows the B 2 /87r contribution to the total pressure. The inset in the lower left panel shows the 
effective mass as a function of B*. In the lower right panel, the electron and neutron concentrations have been 
suppressed for clarity (Y e = Y p — and Y n = 1 — Y p ). These results include the effcts of L anda u quantization 
and the anomalous magnetic moments in interacting nucleon matter. Figure taken from Ref. [159 ] . 



attainable with a magnetic field governed by a constant current function is noticeably larger than that 
attained by rotation. Hydrostatically stable configurations (of which some may not be stable to dynamical 
perturbations) are contained between the heavy solid lines. The lower heavy solid line in each panel is 
the usual field-free, spherical result for the mass-radius relation. The upper heavy solid line represents 
the largest possible stable mass for a given equatorial radius as the internal magnetic field strength is 
increased. Large axialiy-symmetric fields tend to yield flattened configurations, and if large enough, shift 
the maximum densities off-center (these configurations are contained in the yellow s haded regions in the 
figure). This results in toroidal shapes with low-density centers. As discussed in Ref. 1471 ]. regions to the 
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FIG. 16: Limits to hydrostatic configurations for neutron stars permeated by axially symmetric magnetic fields. 
The upper (lower) panel is for the EOS GM3 including (excluding) hyperons. In each panel, the lower heavy solid 
curve is the standard mass-radius relation for field-free stars. The upper heavy solid curve represents the largest 
gravitational mass possible for a given equatorial radius as the magnetic fields are increased, for the indicated 
current function /. The yellow shaded region contains configurations in which the maximum density point is 
off-ce nter . The heavy dashed curve is the locus of minima of the contours of fixed baryon mass configurations: 
Ref. [1471 ] suggested that this could be the limit to dynamical stability. Thin solid lines are contours of the 
maximum magnetic field strength in units of 10 18 G. Similarly, thin dashed lines are contours of maximum mass- 
energy density in units of the energy density po = p s at the nuclear saturation density. Figure taken from [160 ] . 



left of the heavy dashed line, which are the loci of minima of fixed baryon mass configurations, are likely 
unstable to small amplitude perturbations. 

Superimposed on Fig. [Tniare contours of maximum mass density and maximum magnetic field strength 
within a star. For the GM3 model, hyperons appear at zero field at about twice nuclear density. Therefore, 
the portions of the two panels in which the maximum density is below about 3 times nuclear density are 
nearly identical. In the case of the GM3 EOS, with or without hyperons, the maximum field value for any 
presumably stable configuration is B ~ 1.8 x 10 18 G. For a constant current function, and for a variety of 
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EOSs, Ref. [l47j | found that the maximum value of the magnetic field in stable stars never exceeds 3 x 10 18 
G, which corresponds to B* ~ 7 x 10 4 . This field strength is not nearly large enough to produce appreciable 
effects on the EOS. This includes changes to the hyperon or nucleon compositions, with or without the 
inclusion of anomalous magnetic moments. Additionally, the magnetic fields with the assumed current 
function have relatively small spatial gradients so that the ratio of the maximum field to the average 
field within the star is not large. Therefore, the dominant effect of the field arises through the magnetic 
field stress B 2 /8tt, which effectively dominates the matter pressure below a few times nuclear saturation 
density, depending on the field's orientation. Whether other choices of current functions, or the relaxation 
of the condition of axial symmetry, will alter these conclusions has not been explored. 

It might be that the shapes of the stars will significantly change with a different field geometry. It is 
even possible to imagine a disordered field for which (B 2 ) is significantly larger than (B) 2 . In this case 
the pressure will be dominated by fluctuations in the field, but the stars will tend to be spherical. It 
is possible that strong magnetic fields may be held in the core for periods much longer than the ohmic 
diffusion time due to interactions between the magnetic flux tu bes a nd the vortex tubes expected to be 
present in a superconducting or superfluid rotating neutron star 164| . Although the results to date imply 
that average fields within a neutron star cannot exceed 2 — 3 x 10 18 G before stability is compromised, too 
small to have a significant impact on the EOS, it is possible that these conclusions are not true in some 
magnetic field geometries. This is an important subject that deserves additional investigation. 



X. GRAVITY WAVE EMISSION FROM BINARY MERGERS 

Mergers of compact objects in binary systems, such as a pair of neutron stars (NS-NS), a neutron 
star and a black hole (NS-BH), or two black holes (BH-BH), are expected to be prominent sources of 
gravitational radiation [165] . The gravitational-wave signature of such systems is primarily determined 



by the chirp mass M chirp = (MiM 2 ) 3/5 (Mi + M 2 )~ 1/5 , where M 1 and M 2 are the masses of the coalescing 
objects. The radiation of gravitational waves removes energy which causes the mutual orbits to decay. 
For example, the binary pulsar PSR B1913+16 has a merger timescale of about 250 million years, and 



the pulsar binary PSR J0737-3039 has a merger timescale of about 85 million years [32|, so there is ample 
reason to expect that many such decaying compact binaries exist in the Galaxy. Besides emitting copious 
amo unts of gravitational radiation, binary mergers have been prop osed as a source of the r-process elements 
166 ] and the origin of the shorter-duration gamma ray bursters |l67l . 



From the point of view of neutron star structure, observations of gravity waves from merger events 
offer a unique opportunity to simultaneously measure masses and radii, and could set extreme limits on 
the neutron star maximum mass. Current and detectors may not have the ability to receive information 
expected at high frequencies from NS-NS mergers, so our discussion focuses on BH-NS binaries. From the 
inspiral waves waves, masses of the binary components can be measured. The inspiral terminates either 
due to tidal disruption or to the stars' reaching their mutual ISCO. In either c ase, measurement of the 



gravity- wave frequency at this point can be used to infer the neutron star radius [1681 ] . Furthermore, from 



the ringdown phase, the mass of the final merged remnant can be found. The binding energy thus observed 
can be connected to M/R and/or the EOS (see §VIII Bj) . Moreover, a signal might be observable from a 



super-sized neutron star formed from the pair if a black hole doesn't form. This would set a lower limit to 
the neutron star maximum mass that is substantially larger than at present. Importantly, as demonstrated 
below, the wave pattern emitted if one of the components is a strange quark star is substantially different 
than that of a normal hadronic star. 
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A. Tidal disruption and mass transfer 

A neutron star in a binary merger with a larger mass companion will be tidally disrupted. This occurs 
when the radius of the less massive component exceeds the innermost equipotential surface (i.e., its Roche 
lobe) surrounding the two stars that includes the first, or inner, Lagrange point. If the tidal disruption 
occurs far enough outside the innermost stable circular orbit (ISCO) of th e binary, it is e xpected that an 
accretion disc will form or that mass transfer to its companion will occur (l^.[l70l.[l7l|. In either case, 



the gravity wave signal is expected to be significantly different than if tidal disruption occurs after the 
star penetrates the ISCO. 

Tidal disruption or the onset of mass transfer depends upon the EOS of dense matter and the mass ratio 
q = M\jMi of the binary. If mass can be transferred quickly enough from the lighter to the heavier star, 
conservation of mass and orbital angular momentum can be assumed, which generally leads to a reversal 
of inspiral. The binary separation widens. Although this increases the Roche lobe volume, the radius of 
the neutron star also increases in response to the mass loss. Mass transfer will continue in a stable fashion 
if the lighter star can expand sufficiently fast such that it is able to continuously fill its Roche surface. 
Mass transfer under such conditions is termed stable mass transfer. Because the orbital separation now 
increases, the gravity wave amplitude and frequency will decrease. The signature of stable mass transfer 
in gravity waves should therefore be strikingly different than for an amorphous tidal disruption or a direct 
plunge. 

As we show below, important information about the radius of a neutron star and the underlying EOS 
could be contained in the gravity wave signal. For example, the expected gravity-wave signal from a 
merger involving a neutron star should be quite different than from a self-bound strange quark matter 
star. In addition, the longer timescale produced by stable mass transfer might also extend the duration of 
the event from millisecond timescales, the orbital timescale near the ISCO, to a few seconds, which then 



could explain the duration of short-timescale gamma ray burst |171 



Another effect of stable mass transfer would be to modify the amount of material potentially ejected 
from the system. Matter from a tidally disrupted neutron star, which could be accelerated to escape 
velocities from the binary |l66| |. undergoes decompression which re sult s in h eavy nuclei and an intense 
neutron flux leading to the copious production of r-process elements In the case in which stable 



mass transfer occurs, sudden disruption of the neutron star near the last stable orbit is avoided, but mass 
could be ejected more easil y at later times and at larger separations when the neutron star approaches its 



minimum stable mass 1174 



B. Merger evolution 



Here we explore the qualitative character of a binary merger using a relatively simple model that ignores 
relativistic corrections to the Roche ge omet ry and the orbital evolution but still reproduces the qualit ative 



features of more sophisticated models [175 



This model is similar in some respects to that of Ref. |l7ll 
but correctly includes the effects of the EOS. To be definite, we consider a low-mass BH-NS system, with 
M^s — 1.4 M Q and Mbh — 3 — 1OM . The original orbit decays by the emission of gravitational radiation. 
Denoting the total orbital mass M + Mbh + M^s, the reduced mass \i = Mbj jMns/M, and the mass 
ratio q = M^s/Mbh, the rate of change of orbital angular momentum J is |l76 



Jew 



32G 7 /Vikf 5 / 2 32 G 7/2 q 2 M 9 / 2 



5 c 5 a 7 /2 
where a is the orbital semi-major axis and 

J 2 = GMfi 2 a = GM 3 aq 2 (l + q) 



5 c 5 (l + g) 4 a 7 / 2 



119) 



(120) 
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The orbital frequency is 



IGM 



'121) 



These equations are valid for circular orbits, but it can be shown that the timescale for decay of orbital 
eccentricity is much shorter than the timescale for decay of semi-major axis a, which for circular orbits is 
equivalent to the orbital separation. 

The binary shrinks until tidal disruption ensues. This occurs when the neutron star fills its Roche 
radius, the gravitational equipotential surface that passes through the inner Lagrange point, and mass 
begins to flow t o its companion. In Newtonian gravity, the Roche radius Re is well approximated by 
Kopal's formula [1771 

R e /a = 0A6[q/(l + q)] 1 / 3 , (122) 



or a better fit by Eggleton [178 



Re/a = 0.49/[0.6 + q~ 2/3 ln(l + q 1 ' 6 )} . 
Thus, mass overflow begins at the moment that Re = R, or 

a = RiRe/a)' 1 , 



(123) 



(124) 



where R is the neutron star radius. Note that both Re /a and a are decreasing functions of q for fixed 
M and J. Mass transfer will continue in a stable fashion if the star's radius, after an increment of mass 
transfer, is less than Re, so that continued gravitational radiation will result in renewed mass transfer, or 
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where M is assumed constant. Eq. (j!25p defines the parameter a, which is a function of the EOS and M 2 . 
It is a which determines both the onset of stable mass transfer and the subsequent evolution. If Kopal's 
formula is used, the second term on the right-hand side of Eq. (11251) is just 1/3. 

If the mass transfer is assumed to conserve orbital angular momentum, the evolution of the system 
will be defined by 

din J 1 cilna , 1-qdlnq j GW 32 G 3 qM 3 



dt 2 dt 1 
Combining this with Eq. (j!25p . one finds 
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Since both Jew and q are negative, the condition for stable mass transfer becomes 

.dlniRe/a) 
a + 2 (1-q - 1 + g \ e/ ' >0 

ding 
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Using Kopal's formula, Eq. (I122p . this is simply a > 2q — 5/3. For hadronic stars near 1.4 M , a ~ 0, 
so q < 5/6 is the condition. For strange quark stars not near the maximum mass, a ~ 1/3, and q < 1 
is the condition. Thus, all NS-NS and BH-NS binaries containing strange quark matter stars will likely 
have epochs of stable mass transfer. However, some some NS-NS binaries with q ~ 1 may not show this 
phenomenon. 
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During stable mass transfer, a > and, hence, the binary spirals apart. However, if the condition Eq. 
(11281) is subsequently violated, mass transfer becomes unstable at that point and the remaining neutron 
star will quickly tidally disrupt. 

Equations (11251) (with the equal sign) and (11261) determine the time evolution of a and q. From these 
quantities, all other relevant quantities can be determined during mass transfer. The details will depend 
upon the EOS. The two main observables will be the gravity wave amplitude and the frequency of the 
gravity waves as a function of time. The scalar gravitational polarization amplitude (or dimensionless 
strain) is 

4 AG 2 M 2 a 

h ^=r * ( 1+ V °^('- r >. < 129 > 

where r is the distance from the binary system to the observer. The observed frequency of the emitted 
gravitational waves will be twice the orbital frequency u Eq. (I12ip . Typical values are \h + r\ ~ 2 x 10~ 19 
pc and v = uj/2\pi ~ 0.1 — 0.5 kHz. 



C. The EOSs of normal versus self-bound stars 

It is useful to examine the behavior of the function a(M^s for both neutron stars and self-bound 
trange quark matter stars. For the discussion at hand, we will use the term normal star to refer to a star 
with a surface of normal matter in which the pressure vanishes at vanishing baryon density. The interior of 
the star, however, may contain any or a combination of exotica such as (1) strangeness-bearing matter in 
the form of hyperons, kaons, or quarks, (2) Bose (pion or kaon) condensed matter, and (3) quark matter. 

A self-bound star, as exemplified by Wit ten's conjecture that strange quark matter might be more 
stable than normal matter at zero pressure [H, 179 ], has a bare quark matter surface in which the pressure 



vanishes at a finite but supra-nuclear baryon density. In the context of the MIT bag model with first order 
corrections due to gluon exhange, the baryon density at which pressure vanishes is 

n b (P = 0) = (4B/3vr 2 / 3 ) 3 / 4 (l - 2a c /vr) 1 / 4 , (130) 

where B is the bag constant and a c = g 2 /(4:7i) is the quark-gluon coupling constant. This density is not 
signi ficantly affected by the finite strange quark mass 11801 or by the pairing phenomenon in quark matter 



[Tii 



Examples of radius versus mass for broadly differing EOSs selected from are shown in Fig. [171 
While q uant itative differences exist among normal stars, a rough representation of o:(Mns) for low mass 



stars is 1175 



0.09 

a ~ ~ M NS /M Q -0.09 ' ^ ' 

This relation is useful in understanding detailed numerical simulations. Qualitative differences in the 
outcomes of mergers with a black hole emerge, however, because of the gross differences in the mass-radius 
diagram of normal and self-bound stars. Thus, a normal star and a self-bound star represent two quite 
different possibilities (see the right panel in Fig. [1] 



_ dYnR J < for a normal neutron star (NS) (-\oo\ 
a = d\nM 1 > for a self-bound SQM star ' ^ ' 



We will explore the astrophysical consequences of these distinctive behaviors in BH-NS mergers. 
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FIG. 17: Radius versus mass (left panel) and its logarithmic derivative (right panel) for prototype EOSs. The 
EOS symbols are as in [fj]. The vertical lines in the left panel and the horizontal lines in the right panel are 
discussed in text. In the left panel, the maximum masses for which stable mass transfer is allowed in the case 
q = Mns/Mbh = 0.5 are denoted by filled circles for each EOS. 

Note that a is intimately connected with the dense matter EOS, since there exists a one-to-one corre- 
spondence between R(M) and P(ns), where P is the pressure and ub is the baryon density. Gravitational 
mergers in which a compact star loses significant mass during evolution is one of the rare examples in 
which the R versus M (or equivalently, P versus n#) relationship of the same star is sampled. Although 
we focus here on the coalescence of a compact star with a BH, the theoretical formalism and our principal 
findings apply also to NS-NS mergers. 



D. Model evolutions 

Figure [18] shows the results of integrating Eqs. (I125P and (I126P for simulations involving a 4.5 M Q 
black hole and a 1.5 M normal or self-bound quark star. Although both simulations result in stable mass 
transfer, there is a pronounced qualitative difference between the results after mass transfer begins. Inspiral 
is characterized by increases in the orbital frequency u> and scalar gravitational polarization amplitude h + , 
a decrease in orbital separation a, and a fixed q. Stable mass transfer ensues at the "kinks" visible in the 
evolution of these quantities. Within each class of EOS, i.e., normal and self-bound stars, variations in the 
EOS only qualitatively alter the results. During stable mass transfer, the decrease in orbital separation 
and rise in frequency and waveform amplitudes are reversed. 

In the neutron star case, the star spirals outwards and loses mass at an approximately constant rate. 
Mass loss continues for a few seconds until the star approaches but does not decrease below its minimum 
mass (about 0.1 M ). Mass transfer becomes unstable at this point, and full tidal disruption including 



explosive decompression of the remnant is to be expected in this case [174]]. On the other hand, the quark 
star loses mass exponentially and remains at a nearly fixed orbital separation during this period. Mass loss 
continues for a virtually infinite time. As a result, the temporal behavior of the gravity wave emissions, 
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FIG. 18: Schematic behavior of physical and observational variables in mergers between low-mass black holes and 
neutron stars or self-bound quark stars. The total system mass is 6 M and the initial mass ratio is q = 1/3 in 
both cases. The initial radii of the neutron star and quark star were assumed to be equal. The time scales have 
arbitrary zero points. Upper panel displays semi-major axis a (thick lines) and component mass M/vs, Mqs (thin 
lines) evolution. Lower panel displays orbital frequency v (thick lines) and strain amplitude \h+r\ evolution. In 
both panels, solid curves refer to the neutron star simulation and dashed curves to the quark star simulations. 



both of the amplitude and the frequency, from these two types of mergers significantly differ, as seen in 
Fig. [M 

The major effect of incorporating ge nera l relativistic corrections to the potential is to speed up the 



evolution relative to the Newtonian case [1751 ] . Stable mass transfer thus begins earlier in these cases. GR 
corrections also result in a somewhat larger value for the orbital separations following the onset of mass 
transfer. 

Using the relations established in the previous sections, one can readily understand qualitatively and 
quantitatively the results shown in Fig. [THJ For example, the binary separation when mass transfer begins 
is, from Eq. (I124p . a ~ i?(4/3) 1//3 /0.46 ~ 44 km for the two cases shown (assuming R = 12 km). The 
orbital frequency will be a maximum here, v = (oj/2h) ~ 500 Hz, from Eq. (I12ip . 

In the neutron star case, mass loss continues until the condition Eq. (11281) is violated. This occurs 
before the minimum mass is reached, since near the minimum mass a — > — oo. The binary separation 
increases since, using Kopal's formula, Eq. (11221) . dlna/dlnq = (a — l/3)/(l + q) is negative for a < 1/3, 
which is always true for neutron stars. 

However, in the quark star case, the binary separation remains nearly fixed since a — 1/3 ~ 0, especially 
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for q — > 0. The condition Eq. f 1 1 2 8 j) is never violated for quark stars, so mass transfer is not terminated. 
Using Kopal's formula and setting a = 1/3, one arrives at the evolution equation 



dq 
dt 



32 G 3 M 3 q 2 /0.46\ 4 ( Mq 



5 c 5 1 - q 2 I Rq 



V M 



4/3 



(133) 



where Rq and Mq are the quark star's initial radius and mass, respectively. This is trivially integrated to 
find the elapsed time to go from an initial mass ratio qt to a final mass ratio g/ 



At 



Rq\ 



M 



4/3 



32M 3 G 3 V.46/ \M Q> 



0.002g7 1 s. 



(134) 



where we assumed in the final step that qj « qi. Thus, the mass of the self-bound star dwindles to 
extremely small values < 10~ 30 M Q , the minimum mass limit being that of a strange quark nugget, 



determined in part by surface and Coulomb effects. This takes, from Eq. (I134p . more than 10 27 s. 

The envelope of the gravitational waveform amplitude |^+(t)| follows the behavior of q and a. In the 
normal neutron star case, following the onset of stable mass transfer, the amplitude decreases steadily 
until stable mass transfer terminates, but decreases much more rapidly in the the quark star case. 



E. Observational parameters 



As discussed in Ref . [1821 ] , careful analysis of the gravitational waveform during inspiral yields values 
for not only the chirp mass M chirp = (M BH M NS ) 3 ^ 5 /M 1 / 5 , but for also the reduced mass M BH M N s/M, 
so that both Mbh and M^s can be found. Thus, observation of stable mass transfer effects in the 
gravitational wave signal will allow several details about neutron star structure to be discerned, including 
several features of the underlying EOS. For example, the onset of mass transfer can be determined by 
the peak in u, a nd th e value of u there gives a. The Roche limit condition Eq. (11241) . as modified by 
general relativity |l75| . the n all ows the determination of the star's radius. Thus a point on the mass-radius 



diagram can be estimated [183 



This information can be supplemented using the gravitational wave amplitudes h + at its peak and at 
the end of stable mass transfer. The combination h + uj~ l l 3 depends only on a function of q, so the ratio 
of that combination and knowledge of should allow determination of qf. From the Roche condition and 
knowledge of a/ from Uf, another mass-radius combination can be found. 

The time elapsed between the onset of stable mass transfer and its termination depends in a compli- 
cated way on q, M and a(M). Nevertheless, measurement of this time, assuming that q and M are already 
measured, implies that the function a can be constrained. 

Most importantly, the sharp contrast between the evolutions during stable mass transfer of a normal 
neutron star and a strange quark star should make these cases distinguishable even if other information 
concerning the values of q, M and R is weak. This result is independent of the form of the gravitational 
potential or the nuclear or quark matter EOS employed. 

Finally, for the case of strange quark matter stars, the differences in the height of the frequency peak 
and the plateau in the frequency values at later times are related to the differences in radii of the stars at 
these two epochs. (In Fig. [T5J the peak and the plateau are coincident because a = 1/3 was assumed for 
all masses.) It could be an indirect indicator of the maximum mass of the star: the closer is the star's mass 
before mass transfer to the maximum mass, the greater is the difference between these frequency values, 
because the radius change will be larger. Together with radius information, the value of the maximum 
mass remains the most important unknown that could reveal the true equation of state at high densities. 
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XI. CONSTRAINTS FROM LABORATORY DATA 

We turn now to opportunities afforded by laboratory experiments to constrain the dense matter equa- 
tion of state. Supplements to our brief account here are contained in several reviews, some of which are 
alluded to below. 

A. Nuclear masses 

The most fundamental property of a nucleus is its binding energy, or mass defect relative to the total 
individual masses of its Z protons and TV neutrons. The first giant leap to understand the nuclear mass 
systematics was taken by Bethe and von-Weizacker who developed the semi-empirical mass formula using 
the liquid-drop model. Incorporatin g rea listic nuclear surfaces and effects stemming from neutron-proton 



asymmetries, Myers and Swiatecki [185| formulated the droplet model approach in which the nuclear 
energy can be written as 



2 



z 2 z 4/3 

+E dif — + E ex ^ + aAA~ 1/2 . (135) 

In this expression, B ~ 16 MeV is the binding energy per particle of bulk isospin symmetric matter at 
saturation, E s , Ec, E^f and E ex are coefficients for the surface energy of symmetric matter, the Coulomb 
energy of a uniformly charged sphere, the diffuseness correction and the exchange correction to the Coulomb 
energy, respectively. The last term represents pairing corrections, where A is a constant and a = +1 
for odd-odd nuclei, for odd-even nuclei, and —1 for even-even nuclei. For simplicity, the effects of 
curvature and higher-order terms are neglected. The quantity S v is the usual symmetry energy coefficient 
and the quantity S* is related to the surface tension associated with the asymmetry parameter 5 = 
{n n — n p ) I {n n + n p ) according to 

/ 3 \ 2 / 3 

S; = 4tt — - a s . (136) 



Airn s 

Note that the surface energy E s = 47rrQcr(<5 = 0) and the surface symmetry energy S s = Aureus, where 
(47rrQ/3)(0.16 fm~ 3 ) = 1, so that S* and S s differ only in the values of the equilibrium densities used. 

For the description of nuclei in the equation of state relevant for astrophysical simulations of supernovae 
and neutron stars, the droplet approach for isolated nuclei has been extended to the case in which nuclei 
are i mmersed in a dense medium comprised of electrons, positrons, protons, neutrons and alpha particles 



186[. The determination of the various parameters entering Eq. (11351) is afforded by the experimental 
nuclear masses. Fits of Eq. f II 3 5 j) to currently available masses have revealed interesting linear correlations, 
particularly for S s /S v with S v . Theoretically, S s /S v is closely connected to the neutron skin thickness 



5R = (r^) 1 ^ 2 — (r 2 ) 1 ^ 2 . (For recent accounts, see |l87l . 188l | and references therein.) For example, and 



depending on the precise way in which the nuclear surface and Coulomb attributes are treated, correlations 
emerging from equivalent fits to nuclear masses are 

SJS V = -5.253 + 0.254^ or S s /S v = -3.453 + 0.163^ . (137) 

Tighter constraints will be provided by a larger number of experimental masses for nuclei with N > Z, 
highlighting the importance of new mass measurements using rare isotope accelerators (RIA). Fo r new 



developments in techniques of measuring masses of short-lived nuclei off the stability line, see Refs. [189 
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B. Nuclear matter compression modulus 



The task of determining the isospin symmetric nuclear matter compression modulus 

tf(E/A) 
• drfl 

where E/A is the energy particle (determined from the empirical nuclear masses) at the nuclear equi- 
librium density n s ~ 0.16 MeV fm -3 (determined from electron and hadron scattering experiments on 
nuclei), has been arduous on both experimental and theoretical fronts. Thanks to vastly improved exper- 
imental techniques, the identification and analysis of isoscalar giant monopole and isoscalar giant dipole 
resonances have been accomplished for light to heavy nuclei such as 70 Zr, 116 Sn and 208 Pb. Equally chal- 
lenging theoretical analyses of the data based on fully self-consistent Hartree-Fock plus random phase 
approximation (RPA) calculations have now been performed utilizing both non-relativistic potential and 
relativistic field-theoretical models. The conclusions that emerge from these combined experimental and 
theoretical efforts are: 



(138) 

n D 



1. The isoscalar giant monopole resonance data yields K = 240 ± 20 MeV. The uncertainty of about 
20 MeV is mainly due to the density dependence of the symmetry energy around n s ; and 

2. The isoscalar giant dipole resonance data tend to point toward lower values of K. However, there is 
consensus that the extraction of K in this case is largely experimental as the maximum cross section 
decreases very strongly at high excitation energy and likely drops below the current experimental 
sensitivity for excitation energies above 30 and 26 MeV for 116 Sn and 208 Pb, respectively. With 
improved experimental techniques and analysis, these difficulties may be overcome. 

An excellent account o f the se developments can be found in the recent short review (and in other 
articles) by Shlomo et. al. 19Cll . 



C. Symmetry energy and giant dipole resonances 



The analyses of giant resonances, particularly dipole reson ance s, have long served to delineate the role 
of volume and surface effects in nuclei; for a review see Ref. |l91 |. In medium to heavy nuclei in which 



magnetic contributions are small, the inverse-energy-squared weighted photoabsorption cross section cr_ 2 
can be related to the static polarizability p as 

a„ 2 = / du = 27T 2 (e 2 /hc)p . (139) 
J cu 2 



In nuclei with mass number A > 100, er_ 2 = (2.9 ± 0.2) A 5 / 3 /ib MeV" 1 parameterizes the data [l92l. fl93l 



The dipole polarizability can be evaluated as the response of a nucleus to an external dipole field t]D, 
where rj denotes the strength and D = (1/2) J2f=i z i T f is the dipole operator. Explicitly, 

P = 2]T (140) 

where |n) and uo n are the eigenstates and eigenenergies of the nuclear Hamiltonian responsive to the dipole 
operator D. 
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Microscopic RPA calculations [1941 ] of p employing Skyrme-like interactions reproduce the general 
trends of the data. The interplay between volume and surface effects are, however, difficult to extract 
from RPA calculations. Semiclassical methods, in which the energy density formalism is coupled with 
a hydrodynamic approach to describe collective excitations, have thus been employed to explore how 
the v olum e and sur face symmetry energies, S v q and S s , determine the value of p across the periodic 
table (laj Eal, H]. start ing from a given Hamiltonian density, and writing the transition density as 
8[Pn(r) — Pp{ r )] — V ( f ) { r ) cos 6*, the polarizability is calculated from 



p — — / 4>(r)r 3 dr 
3 



(141) 



where 0(r) is the solution of the corresponding Euler-Lagrange (integro-differential) equation. Investiga- 
tions have found that similar results for p are obtained for correlated values of S v and S s . Numerical 
results for 40 Ca, 120 Sn, and 208 Pb show that the choice of S v in th e ran ge 27-42 MeV requires \S S /S V \ = 



1.2-2.2; lower values of S v demand lower values of S s for good fits [196 



A qualitative understanding of this correlation can be gained by using a schematic symmetry energy 



density functional and a leptodermous expansion of the density [191|, whence one obtains the result 



P 



A (r 2 ) 



(142) 



where (r 2 ) is the mean square radius of the nucleus and higher order terms contain corrections from the 
diffuseness and t he skin thickness. Although values of S v = 32.5 MeV and \S S /S V \ = 2.2 describe the data 



adequately [191 
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correlated variations in these numbers are allowed as pointed out in Ref. 
Sum rules have been particularly useful to relate experiments and t heor y in discussing the mean 
excitation energies, widths, and the spreading of the excitation strengths 
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The moments m p of the 



strength function S(u) = J2 n >o \ ( n \F\ 0)1 ~ u n) defined by 



POO 

/ S{u)u p du = J2 \( n \ F \°)\ 

"' n>0 



0J* 



(143) 



where F is the physical operator exciting the nucleus from its ground state |0) to its eigenstate \n), are 
especially helpful in this regard. For example, a good measure of the mean excitation energy is provided by 
E(D) 
data I191H 
Ref. 



'mi/m_i, for which results from RPA and hydrodynamic calculations have been compared with 
Using a droplet model coupled with a hydrodynamic approach to excite the dipole resonance, 



197 obtains 



E(D) 



6h 2 (l + K 



D 



S v 



\ M(r 2 ) [1 + 5S./(S V AW)) 



(144) 



where M is the nucleon mass. The quantity Kd is a model dependent enhancement factor characterizing 
the relative contribution of the nuclear interaction to the mi sum rule and depends critically on the value 
of the energy up to which the energy integration is carried out in analyzing the experimental data. Using 
values of S v = 32.5 MeV and S s /S v ~ 2 fro m a droplet model fit to nuclear energies, and K D = 0.2 
corresponding to E max = 30 MeV, Ref. [l97| finds that a hydrodynamic approach is able to reproduce 
the experimental mean excitation energy of nuclei ranging from 40 Ca to 208 Pb reasonably well. Several 
Skyrme interactions with different values of S v and S s /S v are also able to account for the data to the same 
level of accuracy. However, since the values of S v and S s for these forces are correlated because they are 
fit to experimental masses, the additional constraint on the permitted ranges of S v and S s resulting from 
fitting dipole resonances does not seem to be significant. 
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It must be emphasized that while a hydrodynamic approach is able to account for the gross features 
of nuclear mass dependence, several detailed features of the data such as strength fractionation, spreading 
widths, etc., are not naturally incorporated in its scope. For such details one must adopt a more microscopic 
approach that include s, fo r example, contributions from 2p-2h excitations, etc. For a detailed account, 
see, for example, Ref. |l98 . 



D. Neutron skin thickness in nuclei 



Unlike proton distributions, neutron distributions in nuclei have remained uncertain to this date. 
Studies of neutron densities from a global analysis of medium-energy proton scattering on 208 Pb indicate 
that 0.07 < 5R < 0.16 fm " 
atom data that gives 5R = 



199 . 



Related in form ation is also available from an analysis of antiprotonic 
0.15 ± 0.02 fm |200j ]. In the latter work, nucleon density distributions are 



parameterized by Fermi functions and it is found that the half-density radii for neutrons and protons 
in heavy nuclei are the same, but the diffuseness parameter for the neutrons i s lar ger than that for the 
protons. Skin thicknesses as large as 0.2 fm were obtained in earlier analyses [201 ] . Since these studies 



involve strongly interacting probes, even to this date the value of 5R for a nucleus such as 208 Pb is not 
accurately known. This situation should improve as it is expected that the neutron rms radius will be 
determined to about 1% accuracy by measuri ng th e parity- violating electroweak asymmetry in the elastic 



scattering of polarize d ele ctrons from 208 Pb 2021 



This experiment, named PREX, is planned at the 
Jefferson Laboratory [203l ] for the summer of 2008. 

The neutron skin thickness is of interest as the pressure of neutron star matter below and above the 
saturation density n s depends o n the den sity dependence of the isospin asymmetric part of the nuclear 
interaction. Typel and Brown |205| have noted that model calculations of the difference between 
neutron and proton rms radii 5R = (r 2 ) 1 ^ 2 — (r 2 ) 1//2 are linearly correlated with the pressure of pure 
neutron matter at a density below n s characteristic of the mean density in the nuclear surface (e.g., 
0.1 fm~ 3 ). The density dependence of the symmetry energy controls 5R (we will call this the neutron 
skin thickness) in a heavy, neutron-rich nucleus. Explicitly, 5R is proportional to a specific average of 



[1 — S v / E sym (n)} in the nuclear surf ace, see Refs. 0,@- 



" / ** & - \ /J ■ I ' L ' J 

Horowitz and Piekarewicz 2061 ] have pointed out that models that yield smaller n eutron sk ins in 



2021 . B . have 



heavy nuclei tend to yield smaller neutron star radii. These authors, along with others 
also pointed out the need for an accurate measurement of the neutron skin. 

For the connection b etwe en isospin asymmetry in nuclei and properties of neutron stars, see the recent 
review by Steiner et al. [1871 . 



E. Heavy-Ion Collisions 



1. Collective Flow 



Nuclear collisions in the range E^/A = 0.5 — 2 GeV offer the possibility of pinning down the equation 
of state of matter above norm al nuclear density (up to ~ 2 to 3n s ) from a study of matter, momentum, and 



energy flow of nucleons [2071 ] . The observables confronted with th eoret ical analyses include (i) the mean 



transvers e m omentum per nucleon (p x )/A versu s rapidity y/y pro j [2081 , (ii ) flow angle from a sphericity 



analysis [209 



(iii) azimuthal distributions |210j |. and (iv) radial flow [21l| . Flow data gathered to date 
are largely for protons (as detection of neutrons is more difficult) and for collisions of laboratory nuclei in 
which the isospin asymmetry is not large. Theoretical calculations have generally been performed using 
Boltzmann-type kinetic equations. One such equation for the time evolution of the phase space distribution 
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function f(r,p,t) of a nucleon that incorporates bot h the mean field U and a collision term with Pauli 
blocking of final states is (see, for example, Ref. 212| ) 



^+V p U-V r f-V r U-V p f 



j3 j3 jo d&NN 

a P2 a P2> du- 



v 12 (2n) 3 5 3 (p + p 2 - Pv - P2>) 



(2vr) 6 J " r *~ r * — rffi 
x [// a (l - - /*) - /i'Ml - /)(1 - / 2 )] • 



(145) 



Above, da NN /dQ is the differential nucleon-nucleon cross-section and t> 12 is the relative velocity. In general, 
the mean field U depends on both the density n and the momentum p of the nucleon. Equation f)145|) 
contains effects due to both hard collisions and soft interactions, albeit at a semiclassical level. Theoretical 
studies that confronted data have thus far used isospin averaged nucleon-nucleon cross sections and mean 
fields of symmetric nuclear matter. It is now well established that much of the collective behavior observed 
in ex periments stems from momentum dependent forces at play during the early stages of the collision 
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The conclusion that has emerged from several studies is that as long as momentum dependent forces 
are employed in models that analyze the data, a symmetric matter com pression m odulus of ~ 220 — 240 



214j , fits the heavy- ion data 



MeV, as suggested by the analysis of the giant monopole resonance data jl90l , 
as well (see the recent review by Danielewicz et. al. |215| ). 

The prospects of rare isotope accelerators (RIA's) that can collide highly neutron- rich nucl e i ha s 
spurred further work to study a system of neutrons and protons at high neutron excess 
Generalizing Eq. (11451) to a mixture, the kinetic equation for neutrons is 



2161. [2171. [218 



dt 



+ V P U ■ V r f n - V r U 
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i=n,p 



(146) 



where J n describes collisions of a neutron with all other neutrons and protons. A similar equation can be 
written down for protons with appropriate modifications. On the left hand side of each coupled equation 
the mean field U = U(n n ,n p ;p) depends explicitly on the neutron-proton asymmetry. The connection 
to the symmetry energy arises from the fact that U is obtained from a functional dif ferentiation of the 
Hamiltonian density. Examples of such mean fields may be found in Refs. 1331 . 12171 . 12181 ] . Observables that 
are expected to shed light on the influence of isospin asymmetry include neutron-proton differential flow 
and the ratio of free neutron to proton multiplicity as a function of transverse momentum at midrapidity. 
Experimental investigations of these signatures await the development of RIA's at GeV energies. In this 
connection, it will be important to detect neutrons in addition to protons. 



2. Multi-fragmentation 

The breakup of excited nuclei into several smaller fragments during an intermediate-energy heavy- 
ion collision probes the phase diagram of nucleonic matter at sub-saturation density and moderate (~ 
10 — 20 MeV) temperatures. In this region of the phase diagram the system is mechanically unstable 
if (dP/dn) TjX < 0, and/or chemi cally unstable if (dp, p /dx)T,p < 0. (A pedagogical account of such 
instabilities can be found in Ref. [219(1 ) . These instabilities, which are directly related to the symmetry 
energy at sub-saturation densities [2201 ]. are believed to trigger the onset of multifragmentation. Because 
of the instabilities, matter separates i nto c oexisting liquid and gas phases, which each have different proton 



fractions, i.e. "isospin fractionation" [22 lj . This fractionation is observed in the isotopic yields which can 



potentially reveal information about the symmetry energy. Also, the scaling behavior of ratios of isotop e 



yields measured in separate nuclear reactions, "isoscaling" , is sensitive to the symmetry energy [2221 . 1223 



This scaling is expressed in the empirically observed ratio of fragment yields from two similar systems 
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with different neutron-to-proton ratios: 

Y 2 (N 1 Z)/Y l (N, Z) ~ exp aN+l3Z , (147) 

where the constants a and (5 can be related to the neutron and proton chemical potentials in a canonical 
ensemble description, and thus to the proton fraction of the s ource . To date, there are many suggestions 



of how the symmetry energy may affect multifragmentation [224] . Ongoing research is concerned with 
an extraction of reliable constraints on the symmetry energy from the presently available experimental 
information. 



3. Isospin diffusion 



Isospin diffusion is the process in which the symmetry energy drives the exchange of neutrons and 
protons between nuclei in a heavy-ion collision (see Fig. [T9l . This diffusion process tends to force the 
isospin asymmetry, i.e. the charge to baryon ratio, of the post-collision target and projectile nuclei to be 
equal. The i sosp in content of the post-collision projectile-like fragment have been recently measured in 
experiments |225j |. Starting from a parameterization of the sy mme try energy near the nu clear equilibrium 
density as in Eq. (j4*4"t information from multifragmentation 226 ] and isospin diffusion |227| observables 
have been combined to yield the constraint 0.64 < i < 1.05 and values of 31-33 MeV for S v q, the symmetry 
energy at saturation densi ty. T his constraint is consistent with the Monte Carlo evaluation of the equation 
of state in Akmal et al. [228 ] . computed from two- and three-body interactions which are matched to 
nucleon scattering phase shifts and the energy levels of light nuclei. 
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FIG. 19: Isospin diffusion process during a heavy-ion collision . (Figure courtesy of A. W. Steiner. 



In the case of isospin diffusion, the constraints obtained on the symmetry energy have also led to 
constraints in nuclear structure physics. This simple picture i s, ho wever, modified by fragment emission 
during and after the collision. This complication is alleviated 2291 ] by considering the ratio 



Rg 



25 



A+B 



6 A+A _ S B+ 



B 



§A+A _ fiB+B 



(148) 



where A and B are nuclei with different isospin asymmetries and 5 is the isospin asymmetry of the 
projectile-like fragment. Using experimental data from the collision of 112 Sn and 124 Sn obtained at the 
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National Superconducting C yclot ron Laboratory (MSU) [225| . a value R$ ~ 0.46 was obtained. Isospin- 
dependent transport models |217l | of the same collision suggest that only symmetry energies with 0.64 < 
7 < 1.05 as discussed above are compatible with the observed value of R$. Because the symmetry energy 
is tightly corre lated with the neutron skin thickness (the difference between the neutron and proton RMS 



radii) in lead 204j , the heavy-ion collision data also offer a constraint on the neutron-skin thickness, 



Rp > 0.15 fm 
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Although the connection between the heavy-ion data and neutron star radii is more tenuous, owing to 
the fact that neutron star radii are sensitive to the symmetry energy at higher densities, isospin diffusion 
also offers guidance on the radii of 1.4 solar mass neutron stars. The symmetry energ y from Eq. 
suggests that the radius of a 1.4 solar mass star is likely between 11.5 and 13.6 km 
consistent with presently available astrophysical observations. 
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XII. OUTLOOK 

Even relative to a decade ago, the quantity of data and the variety of approaches using light waves 
of all wavelengths to estimate neutron star masses, radii, compactness, crustal thicknesses, etc. have 
mushroomed. Traditional techniques such as pulsar timing measurements have produced more than two 
dozen masses accurate to the 10% or better level, enough that their statistical analysis are beginning 
to yield information about their origin and evolution. More than a dozen thermally emitting neutron 
stars have had radiation radii estimated. A few stars have had redshifts determined either via direct 
measurements of atmospheric absorption lines or from estimates obtained from light-bending arguments. 
Regions in the M — R plane have been staked out from observations of accreting neutron stars including 
QPOs and Eddington-flux-limited sources. Other observations, including pulsar glitches, seismological 
studies of neutron star oscillations, and the cooling behavior following bursts or accretion, have constrained 
the extent of neutron star crusts. With the exception of masses, none of these other observations are yet 
sufficiently precise or model-independent to significantly impact dense matter theories. 

In fact, some apparently contradictory conclusions have been reported over the last few years. Some 
neutron stars are observed to be spinning rapidly enough to exclude extremely stiff density-dependent 
symmetry energies. But the best-studied cooling, neutron star, the nearby source RX J1856-3754, indicates 
a relatively large radius. Rapidly cooling neutron star crusts following superbursts and the non-detection 
of neutron stars in several young supernova remnants both indicate high rates of neutrino emission in 
their interiors. Yet the majority of thermally-emitting neutron stars are consistent with relatively slow, or 
standard, core neutrino rates. 

There are indications that new observations and techniques will mount even more rapidly than in the 
last decade, which gives hope that some of these contradictory results will be resolved. The number of 
known pulsars is growing rapidly, and new mass measurements implying that the neutron star maximum 
mass is greater than 1.6 M and perhaps as large as 2 M now exist. This is large enough to perhaps 
delimit the amounts of exotic matter (i.e., hyperons, Bose condensates, or deconfined quarks) in neutron 
star cores. Several new pulsars have been found in highly relativistic binaries, including, in some cases, 
other neutron stars or pulsars. The most highly relativistic system discovered to date, PSR 0737-3039, 
might reveal the moment of inertia of one of its stars, yielding perhaps the most accurate measurement 
of a neutron star radius possible. Furthermore, its relative closness and dimness strongly suggests that 
several even more relativistic binaries will soon be found. 

At the opposite extreme, there are indications that neutron stars smaller than 1.2 M Q exist. This 
value is uncomfortably close to the empirical minimum found in supernova models stemming from the 
evolved cores of the smallest stars thought to be capable of the event, approximately 8 M Q initial mass. 
Theoretically, proto-neutron stars as small as 0.9 M could be stable, and it will be interesting to see if 
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masses approaching this value are observed. 

Less traditional methods, involving neutrinos and gravity waves, that could be utilized to infer in- 
formation about neutron stars should soon become available. The 20 neutrinos detected from SN 1987A 
twenty years ago will be dwarfed by the thousands or even tens of thousands that could be observed today 
from a galactic supernova. Statistical uncertainties in the total binding energy, the deleptonization and 
cooling timescales, the average neutrino energies, etc. should be quite small, permitting important conclu- 
sion to be drawn about the proto- neutron star's mass and opacities. In addition, tens or more i/ e 's could 
be detected during the collapse prior to bounce, their number and average energy crucially depending on 
the density dependence of the symmetry energy below nuclear density. The observation of gravity waves 
from compact star binary mergers is eagerly anticipated, even more so now that the pulsar binary PSR 
0737-3039 has been discovered. In fact, the gravity waves observed from a binary merger might be the 
least unambiguous way to learn of the existence of quark matter stars. 

Complementing the increase in observational data has been the development of new techniques, sup- 
porting theory and laboratory data relevant to dense matter physics. It is now known that the underlying 
physics behind the neutron star maximum mass and typical neutron star radii are quite different. The 
maximum mass is set by the dense matter equation of state well beyond 2n s , but neutron star radii are 
determined by dE sym /dn in the range 1 — 2n s . It is similarly expected that better models of emission 
beams from accreting neutron stars, pulsar glitches, superfluidity and superconductivity involving nucle- 
ons and hyperons, neutron star atmospheres, and the seismology of crustal oscillations will make better 
use of accumulating data. 

Laboratory experiments relevant to nuclear astrophysics in general, and to neutron stars in particular, 
have accelerated significantly in recent years. With impressive efficiency, masses of most nuclei between 
the path of beta stability and the neutron drip line have now been measured. With the Rare Isotope 
Accelerator and Jefferson Lab experiments, many more mass measurements of extremely neutron-rich 
nuclei and flow measurements of neutron-rich heavy ion collisions will be forthcoming. More knowledge of 
hyperon-nucleon interactions through studies of hypernuclei should be available. Perhaps one of the more 
intriguing experiments is the PREX experiment at the Jefferson Laboratory to measure the neutron radius 
of Pb 208 to 1% accuracy, which should yield a neutron skin thickness estimate sufficiently precise to offer 
significant constraints on the neutron star radius. Such experiments should be extended to other neutron- 
rich nuclei. Multifragmentation and isospin diffusion studies from heavy ion collisions should yield more 
accurate information concerning the symmetry energy of nucleonic matter and the critical density for phase 
separation in dense matter, which is particulary relevant to the determination of the core-crust interface 
density, pressure and enthalpy. New studies of giant monopole, dipole and other nuclear vibrational 
modes should be undertaken to study the nuclear incompressibility and symmetry energy properties. In 
addition, new neutrino laboratories are planned, including UNO, that will refine measurements of neutrino 
oscillations which could be relevant for interpreting neutrino signals from supernovae. These new neutrino 
laboratories will also improve by an order of magnitude or more, the number of neutrinos observed from 
a supernova, and will increase the rate at which supernovae are observed as well. New experiments 
involving neutrino interactions with heavy nuclei could shed light on neutrino absorption and scattering 
rates, quantifying the effects of density, spin and isospin correlations for opacities of dense matter. 

Although we've not commented on recent developments in microscopic theory, a topic to which Hans 
Bethe devoted much of his early career, there have been many of interest for neutron stars. Much effort has 
gone into determining the EOS of neutron matter at subnuclear densities which is consistent with phase 
shift data. The effects of Pauli blocking and density, spin and isospin correlations on neutrino opacities 
have been explored. New explorations of the pairing gaps of nucleons and hyperons, the hyperon-nucleon 
interaction and the threshold densities of hyperons have appeared. Finally, new studies are targeting 
the properties of pion and kaon condensates and delineating quark matter properties including color 
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compositions, gaps and critical densities. 
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